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Abstract 

Path integral formulations for the Smorodinsky-Winternitz potentials in two- and three- 
dimensional Euclidean space are presented. We mention all coordinate systems which sep- 
arate the Smorodinsky-Winternitz potentials and state the corresponding path integral for- 
mulations. Whereas in many coordinate systems an explicit path integral formulation is not 
possible, we list in all soluble cases the path integral evaluations explicitly in terms of the 
propagators and the spectral expansions into the wave- functions. 
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1. Introduction. 



In the study of the Kepler problem and the harmonic oscillator it turns out that they possess 
properties making them of special interest, for instance, all finite classical trajectories are closed 
and all energy eigen-values are multiply degenerated. These properties point beyond their 
specific shape and give rise to ask about their "hidden" or "dynamical" group structure, namely 
whether it is possible to distort them in a systematic way, and at the same time keep these 
properties. Furthermore, do they belong to a class of potential problems which have similar 
properties? 

For instance, as has been known for a long time, at the basis of the specific properties of 
the hydrogen atom lies its 0(4) symmetry, or to be more precise, the dynamical group 0(4) 
for the discrete spectrum and the Lorentz group 0(3, 1) for the continuous spectrum. All such 
systems have the particular property that all the energy-levels of the system are organized 
in representations of the non-invariance group which contain representations of the dynamical 
subgroup realized in terms of the wave- functions of these energy-levels J53] . In the case of the 



hydrogen it enabled Pauli [112] and Fock [39] to solve the quantum mechanical Kepler problem 
without explicitly solving the Schrodinger equation. Actually, the algebra of the dynamical 
symmetry of the hydrogen atom turns out to be a centerless twisted Kac-Moody algebra j23|] . 

Smorodinsky, Winternitz et al. started a systematic search to find and classify potential 
problems in two and three dimensions which can be seen as non-central generalizations of the 
Coulomb-, the harmonic oscillator and radial barrier potentials. The classification scheme starts 
from the consideration in which way integrals of motion (in classical mechanics), respectively 
additional operators corresponding to these integrals of motion which commute with the Hamil- 
tonian, are related to the separability of the potential problem in more than one coordinate 
system. In two dimensions [42] it turns out that there are four potentials of the sought type 
which all have three constants of motion (including energy), i.e. there are two more operators 
commuting with the Hamiltonian. Smorodinsky, Winternitz et al. extended their investigations 
in a classical paper [93] to three dimensions by listing all potentials which separate in more than 
one coordinate system. However, in three dimensions one must distinguish two cases, namely 
minimally super-integrable systems, where four, and maximally super-integrable systems, where 
five functionally independent integrals of motion exist. 

Actually in D dimensions a system is called minimally super-integrable if it has 2D — 2 
functionally independent integrals of motions, and it is called maximally super-integrable if it has 
2D — 1 functionally independent integrals of motion jj], [34|] . A discussion of the two-dimensional 
case is due to e.g. Winternitz et al. [|129|], and a detailed study of the three-dimensional case can 



be found in [35, 82]. Moreover, the Smorodinsky- Winternitz potentials allow a generalization to 
D dimensions and are then also maximally super-integrable |3q] . 

In the three-dimensional Euclidean space there are eleven coordinate systems which separate 
the free Hamiltonian [74, 101, 105]. Each potential must now be analysed according to its 
separability in these coordinate systems. Because each potential can be rotated and translated 
with respect to itself, one has to look for the equivalence classes among them. As it turns out 
there are eight minimally super-integrable (including a ring-shaped oscillator and the Hartmann 
potential) and five maximally super-integrable (including harmonic oscillator and Coulomb) 
potentials. Evans (35|] presented a list of these potentials including the corresponding integrals 
of motion and (almost) all separating coordinate systems. The group structure of the two- 
dimensional potentials was given in [42], for the Hartmann potential by Gerry [46] and Granovsky 
et al. 5C], for the ring-shaped oscillator by Quesne 116 ] and Zhedanov [131], and the harmonic 
oscillator including l/rr 2 -terms (a Smorodinsky- Winternitz potential, see below) was analysed 
by Evans [35]. 
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The purpose of this paper is to give a path integral discussion of all Smorodinsky-Winternitz 
potentials in two- and three-dimensional Euclidean space. This includes the statement of the 
propagator (the time-evolution kernel), respectively, if only available, the corresponding Green 
function, and the spectral expansion into wave-functions and energy- levels. In the discussion of 
the two-dimensional Smorodinsky-Winternitz potentials the evaluation will be somewhat more 
elaborated in order to demonstrate the techniques. Of particular importance will be the basic 
path integral solutions of the radial harmonic oscillator by Peak and Inomata [ |113| (compare 
also Duru [^] and Goovaerts ^M), and of the Poschl- Teller potential by Bohm and Junker j|] 
which is based on the SU(2) path integral (compare also Duru |28|], and Inomata and Wilson 
|j72|l ). All these path integral solutions are now well-known in the literature and will not be 
repeated here. 

One may ask, if it is possible to analyse the path integral of a potential problem in terms 
of its dynamical symmetry group |fTo[ [71]] . In order to look at such a path integral formulation 
we consider the Lagrangian C(x,x) - V[x) (x € JR p+q ) as formulated, say, in a not-necessarily 
positive definite space with signature 

( gab ) = diag( +1, . , +1 , -!,...,-!) ■ (1) 

p times q times 

One introduces, say, polar coordinates x v = re u (9\, . . . ,6 p+q -i) (y = 1, . . . ,p + q), where the 
e are unit vectors in some suitable chosen (timelike, spacelike or lightlike) set ||. One then 
expresses the Lagrangian in terms of these polar coordinates and seeks for an expansion of the 
quantity e z ^ v ^ 2 > expressed in terms of group elements <?i,<?2- If this is possible one can re- 
express the path integration of the coordinates x into a path integration over group elements g 



yielding |, H @, U4 



x(t")=x" 

- ft' 



x(t')=x' 



J 2?x(t)exp ijf £(x,S)dt i > J dE x d x J2fLDL(9' V) (2) 

= / dE x d x J2fLJ2 D m :k(9')Di <k (g") . (3) 



Here j„ m is defined via the Fourier transformation 

f(g) = [ dE X d X J2fmnDmn(9) , fmn = I S X^U^dg , (4) 

J mn JG 

and dg is the invariant group (Haar) measure. / dE x stands for a Lebesque-Stieltjes integral 
to include discrete (/ dE x — > J2\) as wen as continuous representations. / dE x is to be taken 
over the complete set {A} of class one representations. d x denotes (in the compact case) the 
dimension of the representation and we take 

d\ / Dmn(9) D m>*n'{9)dg = ${\X)8 m ,m'$n,n' (5) 



G 

as a definition for d x . <5(A, A') can denote a Kronecker delta, respectively, a 5-function, depending 
whether the quantity A is a discrete or continuous parameter. We have furthermore used the 
group (composition) law 

DLXda ^b) = E D t(9a)Dl m (g b ) . (6) 

k 

Choosing a basis {6} in the relevant Hilbert space fixes the matrix elements D^ nn through 
Dmn = (D x (g)b m , b n ) of the representation D x (g) of the group. In particular the -Dg m are called 
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associated spherical harmonics, and the Dq are the zonal harmonics. These spherical functions 
are eigen-functions of the corresponding Laplace-Beltrami operator on a, say, homogeneous 
space, and the entire Hilbert space is spanned by a complete set of associated spherical functions 
-Dg m (Gelfand et al. 



j, 44] and Vilenkin | 126 |). In |5£| a systematic study was presented of 
path integrals of the free motion in two- and three-dimensional spaces of constant curvature, 
i.e. on homogeneous spaces. 

The Smorodinsky-Winternitz potentials can also be discussed in the context of a group path 
integration, provided the dynamical "hidden" group structure behind the potential problem 
is known, and the wave-functions and the energy-levels follow from the representations of the 
group [f7(], fn]] . However, things are not as simple as for the free motion because we have not 
a problem of quantum motion on a homogeneous space. Whereas in many coordinate system 
realizations we can actually evaluate the path integral, there remain several ones, where a direct 
path integral evaluation is not possible because we obtain quartic, sextic, complicated Mathieu- 
like path integral problems, where no solution is known. Also a consideration along the lines 
of "quasi-exactly solvable path integrals" [89] can not be applied in an obvious way. However, 
sometimes the zero-energy Green's function can be evaluated. As we will see, the technique 
of separation of variables in the path integral |)5| will give us a convenient tool to exploit the 
actually symmetry structure in such a way that we can systematically reduce the number of 
degrees of freedom by using the appropriate path integral identities. 

We will therefore provide the path integral formulations of all Smorodinsky-Winternitz po- 
tentials in two- and three-dimensional Euclidean space. Some can only be stated, others can be 
explicitly solved. We find that each of the potentials can be at least solved in one coordinate 
system. The identities arising from the relations between the various formulations serve as a 
tool to analyse the "hidden" dynamical group structure, however, only in an implicit way. The 
representations of the group in question where a path integral is not explicitly soluble remain 
still unknown. Furthermore, the explicit spectral expansions corresponding to various coordinate 
systems allow to study the interbasis expansions of the problem. An example of particular im- 
portance is the two-dimensional hydrogen atom [24, 95], respectively the ring-shaped oscillator 
pO|| , an d the Hartmann potential |H], |131| . 

The knowledge of soluble examples and their relation with each other, in particular here the 
Smorodinsky-Winternitz potentials, is of great value in many branches of theoretical physics. 

, pll with either an oscillator or an Coulomb inter- 



For instance, the Hartmann system [33 



action term appears in the study of molecule physics (modeling ring-shaped molecules like the 
benzene molecule), in the study of the Aharonov-Bohm effect [y, |8Q] , for path integral studies 
c.f. e.g. (Bernido and Inomata Cheng [16|, Chetouani et al. |ll|, Gerry and Singh [45], Hoang 
and Giang p||], Sokmen [117], Tanikella and Inomata |123 ], and Wiegel |127|), and for a more 
recent discussion of a deformed generalization [|7(| . Evans [ 35[ also mentioned cases in celestical 
mechanics (e.g. an artificial satellite orbiting an oblate planet, dynamics of elliptical galaxies). 
The consideration of the Smorodinsky-Winternitz potentials in spaces of constant (positive and 
negative) curvature models a space-time of closed, respectively open universes, c.f. the Kepler 
problem and the Higgs oscillator problem (Granovsky et al. [51], Higgs [65], Infeld [67], Leemon 
§7|, Otchik and Red'kov |10§], Pogosyan et al. p5[| , Schrodingerpl]] , and Steveson |12^], and 
path integral discussions c.f. Barut et al. g, ||, and Grosche [Ml ) . Path integral discussions 
of Smorodinsky-Winternitz potentials in space of constant positive and negative curvature will 
be given in forthcoming publications. Starting from a particular exactly soluble example one 
can perturb it by a further (vector or scalar) potential which is not destroying separability in 
a conveniently chosen coordinate system |93[| , one can expand about the exact solution of the 
unperturbed one, hence giving a perturbation series in, say, either powers of the coupling pa- 
rameter (e.g. |100fl ), respectively in powers of h. The latter then gives a "pseudo" semi-classical 
expansion a la DeWitt-Morette @-|03]. 
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Our work can be seen as a generalization and extension of the earlier work of Carpio-Bernido 
et al. [lC]-[p^]. They discussed only some of the Smorodinsky-Winternitz potentials and they 



restricted themselves only to a spherical polar coordinate path integral discussion in the 
three-dimensional Euclidean space. Their results are, of course, included in our discussion. 
They also did not discuss a classification scheme of the potentials. In particular, all our path 
integral evaluations in parabolic coordinates will be new path integral solutions, including the 
expansion into discrete and continuous wave-functions (if present). The path integral solutions of 
the two-dimensional Smorodinsky-Winternitz potential V4 and three-dimensional Smorodinsky- 
Winternitz potential V3 are entirely new (in Refs. [|l0]]-[|l4j] special cases of the latter have been 
discussed). We therefore give a combined presentation of henceforth solutions which have been 
scattered around in the literature and we also give all admissable solutions of each potential 
problem. 

Our paper will be organized as follows: In the next Section we will list for completeness the 
coordinate systems in two- and three-dimensional flat space which separate the free Helmholtz, 
respectively Schrodinger equation. Also a separation formula for path integrals will be stated. 

In the third Section we discuss the Smorodinsky-Winternitz potentials in two and three 
dimensions, respectively. This includes the formulation of the path integral in the various 
separating coordinate systems, and if possible via the basic path integrals, its solution in terms 
of the propagator, respectively the Green functions, and its spectral expansion. In particular, we 
will make frequently use of the separation technique in path integrals as developed in |55[j , and 
will not write down explicitly each step in the calculation of the examples of the explicitly soluble 
Smorodinsky-Winternitz potentials. In our path integral solutions we will state the Feynman 
kernels as explicitly as possible, and their spectral expansions according to 



K(x", x'\ T) = J dE x y* x \(x")V x {x")e- iE * T / h (7) 

f ^n(x')^ n (x")e- iE - T / h 
n 

Jdp%(x f )%(x")e- in ^ 2M (8) 

v n J 

whether the spectrum is purely discrete, purely continuous, or discrete and continuous. In these 
cases, where only the (energy-dependent) Green function is available, the spectral expansions 
will be presented similarly. 

The fourth Section is devoted to a summary and a discussion of our results. 

2. Coordinate Systems in Euclidean Space and the Path Integral. 
2.1. Coordinate Systems in Euclidean Space. 

2.1.1. Coordinate Systems in D Dimensions. In D-dimensional Euclidean space we have the 
coordinate system corresponding to the elliptic coordinates 

IIi#(ei-ei) (9) 
ei < pi < e 2 < ■ ■ ■ < e D < po ■ 

This is the most general coordinate system possible, and all other coordinate systems are in- 
cluded as degenerations. In the following we will explicitly state the coordinate systems in 
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two- and three-dimensional flat space and E^ which separate the Helmholtz, respec- 
tively the Schrodinger equation. We rewrite the Hamiltonian, respectively the Laplacian Al = 
g~ l l 2 d a g ab g 1 l 2 d} > in terms of the corresponding momentum operators p a = —\hg~ 1 ^d a g l / 4: = 
—\h[d a + r a /2) (r a = d a \n.y/g, gdet(g a b)), where we choose a specific ordering prescription 
called product ordering |]52| . Here it is assumed that g a b = h ac h c b, g ab = h ac h cb , the h ac being 
vielbeins, such that we can write H = —t^Al = ^jh ac p a pbh cb + AVpf with the quantum 
potential AVpf given by 

AV PF = ^ [g ab T a T b + 2(g ab T b ), b + g°»^] + ^(2h ac h bc , ab - h a \ a h b \ h - h a \ b h b \ a ) . (10) 

The potentials which may be added without spoiling the separability were derived by Eisenhardt 
[32] and Smorodinsky and Winternitz et al. 93]. 

2.1.2. Coordinate Systems in Two Dimensions. 

Cartesian Coordinates. We consider the usual cartesian coordinates (x,y) = x £ IR 2 = E^ 2 \ 
The metric is given by (g a b) = ^2 and for the momentum operators we have p x = —ihd x , 
p y = —ihdy. Therefore: 

%2 a & ( 92 d "\ 1 -a 

The most general potential separable in cartesian coordinates has the form 

V{x) = u{x)+v{y) . (12) 
Polar Coordinates. Next, we consider two-dimensional polar coordinates 



x = p cos 4> ; P > , 

y = psin</> , < 4> < 2ir . 

The metric is given by (g a b) = diag(l,p 2 ), and the momentum operators have the form 



(13) 



hfdi\ h d 

^=iU + 27J' ^ = 7^- (14) 



This gives for the Hamiltonian 

h 2 fi 2 ( d 2 1 d 1 d 2 \ 1 / 1 \ h 2 

" 2M Ae(2) = "2Mlv + ^ + W = 2Mv' + 7 P V ~ ' (15) 

A separable potential has the form 

V(x) = u(p) + \v(<j>) . (16) 

P 1 

Elliptic Coordinates. Third, we consider the coordinate system 



x = d cosh £ cos n , £ > , 

y = <isinh£sin7/ , —it < rj < ir 



(17) 



and d = R/2, where R is the interfocus distance. The metric is (g a b) = <i 2 (sinh 2 £ + sin 2 17) II2 
and we obtain for the momentum operators 

h / d sinh f; cosh £ \ fi ( d sin n cos r\ \ 

P ^ i \<9£ sinh 2 £ + sin 2 n) ' n i\dn sinh 2 £ + sin 2 r/ ) 
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Consequently we have for the Hamiltonian 



2M E(2) ~ 2Md 2 (sinh 2 e + sin 2 77) Ue 2 + 5?7 2 

= ^^(sinh^ + sin 2 !?)- 1 / 2 ^^^)^^^^^-!^ _ (ig) 

Replacing x ^ x = (f(cosh £ cos rj + 1) yields a coordinate system which we call in the following 
elliptic II coordinates |95[| . The elliptic coordinate system is useful in the study of two center 
systems. Furthermore, an additional parameter is introduced, the interfocus distance R. This 
distinguishes it from the other two-dimensional coordinate systems which can be obtained from 
the elliptic one by considering specific limits. In particular R — » reproduces polar coordinates, 
and R — > oo cartesian coordinates. In the case of the elliptic II systems, R — ► oo reproduces 
parabolic coordinates (see below). 

Note AVpp = which is a peculiarity of two dimensions with metric oc 1, c.f. ([n]). A 
potential separable in these coordinates reads 

= u( C o S H) + v(co SV ) 
sinh £ + sin rj 

Parabolic Coordinates. Finally, we consider the coordinate system 

x = W-f) , y = Zv , £eiR,77>0 (21) 

(alternatively £ > 0, 77 G IR |f)7|| ). We have (<7 a b) = (£ 2 + r / 2 )i2> and consequently for the 
momentum operators 

This gives for the Hamiltonian 

These coordinates have been used in the literature to discuss the two-dimensional "Coulomb- 
problem", c.f. J3l], 42, 57, |68|. The transformation (x, y) 1— > (£, rj) actually is the two-dimensional 
realization of the Levi-Civita transformation. A separable potential has the form 

_ m±^_ (24) 



2.1.3. Coordinate Systems in Three Dimensions. 

Cartesian Coordinates. Again, we start with the simplest case, cartesian coordinates (x, y, z) 
= x £ K 3 = E^ 3 \ Then (g a b) = I3, and p = — iKV. This gives for the Hamiltonian 

" Wl^ = "2^ V2 = 2^ < (25) 
The most general potential separable in cartesian coordinates has the form 

V(x) = u(x) + v(y) + w{z) . (26) 
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Circular Polar Coordinates. Next, we consider circular cylinder coordinates which are very 
similar to the two-dimensional polar coordinates 



x = p cos < 
V = psmc, 
z = z , 



P>0 , 

< <j) < 2tt , 

zeJR . 



(27) 



The metric reads (g a b) = diag(l, p , 1), and the momentum operators are given by (14) and 
p z = —\fid z . This gives for the Hamiltonian 



2M" 



■A 



h 2 ( d 2 Id Id 2 d 2 \ 

EVi) ~~^{d P ^ + 'p'dp' + 'p^'cW + d^ 2 ) 



i 

2M 



Ti 1 



8Mp 2 



A separable potential has the form 



1 



V(x) = u(p) + —zv{4>) + w(z) . 
p 2 

Circular Elliptic coordinates. We consider the coordinate system 



• (28) 



(29) 



x = d cosh £ cos r\ 
y = d sinh £ sin 77 

z = z , 



—ir < r] < it 

zelR 



(30) 



(alternatively £ € IR, < 77 < 7r [|97|]). The metric is (gab) = diag[d 2 (sinh 2 £ + sin 2 rj) , d 2 (sinh 2 £ + 
sin 2 77), 1], and we obtain for the momentum operators dlSj), and p z = —\fid z . Consequently for 
the Hamiltonian 



2M Z 



2 r 



~2M 
1 

2M 



1 



+ 



c) 2 



d 2 (sinh 2 £ + sin 2 77) V <9£ 2 0r/ 2 / d 



+ 



a 1 



1 

d 2 



(sinh 2 £ + sin 2 r])~ 1/2 (p 2 +p 2 )(sinh^£ + sin 2 r/)~ 1/2 + p : 



-1/2 



(31) 



Similarly as in the two-dimensional case, circular elliptic II coordinates are described by the 
replacement x >— * x = d(cosh £ cos 77 + 1). A potential separable in elliptic coordinates reads 



n(cosh£) + 7j(cos77) 
sinh £ + sin 77 



(32) 



Circular Parabolic Coordinates. The last example for cylinder coordinates in three dimensions 
are the parabolic cylinder coordinates 



x 



2 

z = z 



W-e), y = Cv 



£ € M,r] > , 

zeR, 



(33) 



which is the obvious generalization of the two-dimensional case. We have (gab) = diag(£ 2 + 
77 2 ,£ 2 + 77 2 , 1), hence for the momentum operators (22) and p z = —ihd z . This gives for the 
Hamiltonian 



_ ji 2 h 2 

2M E(3) ~ 2M 
1 

2M 

A separable potential has the form 

V(x) 



1 



d 2 d 2 
+ 



+ 



£ 2 + 7/ 2 \ <9£ 2 drj 2 

(e+ifr^Pi+P^e+v 



dz 2 _ 
2^-l/2 



u(£) + v(r]) 



+ w(z) 



(34) 



(35) 
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Sphero- Conical Coordinates. We consider the coordinate system in terms of Jacobi 

x = rsn(p, fc)dn(^, k') , r>0 , k 2 + k! = 1 , 

y = rcn(p, fe)cn(z^, k') , —k < p < k 

z = rdn(p, fc)sn(i/, k') , -2/c' < < 2k' . 

The metric tensor g ao in these coordinates has the form 

(dab) = diag[l, r 2 (k 2 cn 2 p + k' 2 cn 2 i>), r 2 (k 2 cn 2 p + fc /2 cn 2 ^)] , 

the momentum operators are p r = —\h(d/dr + 1/V), and 

fi f d /s 2 sn//cnpdnp \ % / d fc /2 sm/cnz/dnf 

^ i \dp k 2 cn 2 p + k' 2 cn 2 u) ' " i k 2 cn 2 p + k' 2 zv?v 



The Hamiltonian has the form 

-—A -^1 

2M £(3) ~ 2M 



d 2 2 d 1 1 / a 2 ^ 

dr 2 r dr r 2 k 2 cn 2 p + k' 2 cn 2 v \dp 2 du 2 



1 / 2 , 1 1 , 2 , 2 A L 



2M \ r ^k 2 cn 2 p + k' 2 cn 2 v ^k 2 cn 2 p + k' 2 cn 2 v 

Here a separable potential must have the form 

, . 1 f (cna) + u>(cn/3) 
V(x) = u{r) + 



r 2 k 2 cn 2 a + k' 2 cn 2 {3 
The algebraic representation of these coordinates has the form 



x = r H r , r > 



y 2 = -r 2 -^ ^ff^ ^ , oi < pi < a 2 < P2 < a 3 



2 2 

2T = r 



(pi 


-oi)(P2 


- «l) 


(03 


- ai)(a 2 


- ai) 


(Pi 


- Ct2)(P2 


- a 2 ) 


(03 


- a2)(ai 


- a 2 ) 


(Pi 


- a 3 )(p2 


- a 3 ) 


(ai 


- 03) («2 


- 0,3) 



Here we have [P(p) = (p — «i)(p — a 2 )(p — 03)] 



, r 2 p 2 - Pi r 2 Pi ~ P2 
(gab) = drag 1 



4 4 P( P2 ) 



The momentum operators are 



h{ d 1 1 ip'(pi)\ h( d i 1 iP'(p 2 



iV#Pi 2 Pl -p 2 4P( Pl )y ' iV5p 2 2p 2 -pi 4P(p 2 ) 

Together with the zweibeins 

2 V P2 - Pi 2 V pi - p 2 



we obtain for the Hamiltonian 
H = 



tf_(cP_ 2d_ i_ A 
2M \ dr 2 r dr r 2 LB 



1 

2M 



1 



Pr + 12 E h pipi P%h^ + AV PF ( Pl 



. i=i 



8Mr 2 



A 



LB 



92 ~ Pi 











-P(P2 



- d 




-P{P2) 



_d_ 

dp2 



(45) 
(46) 



where the AVpf(Pi) are determined by (|lC|). The relation between the coordinates (pi,/0 2 ) and 
(p, v) is established via 



(at - pi) = (at - a 2 )sn 2 (/x, k) , (a 2 - P2) = (a 2 - a 3 )cn 2 (i/, fc') 
Spherical Coordinates. We consider the spherical coordinates 



(47) 



x = r sin f cos < 
y = r sin sin <j 
z = r cos , 



r > , 

< # < vr , 

< 6 < 2vr . 



(48) 



These are the usual three-dimensional polar coordinates. The metric tensor is (g a b) 
diag(l, r 2 ,r 2 sin 2 9), and the momentum operators have the form 



Pr 



hf d 1 
i \dr r 



Pe 



hfd 1 



hd_ 
i 86 



(49) 



For the Hamiltonian we obtain 



2M 



2M 
1 

2M 



d 2 2d 

dr 2 r dr r 2 \ d9 2 

„2 



+ cot I 



d9 



+ 



2 



r 2 sin 2 1 



, 2 1 2 1 

r A r 2 sm z 



8Afr 2 



1 + 



sin 



(50) 



A potential which is separable in spherical coordinates reads 



1 



V{x) = u{r) + -xv[0) + 



1 



r 2 r 2 sin 2 



-w(6) . 



Parabolic Coordinates. We consider the coordinate system 

y = £rj sin 6 , 



x = £r] cos < 



z 



0< (b<2ir 



This gives for the metric tensor (g a b) 
operators we get 



diag(£ 2 + i] 2 ,( 2 + rj 2 ,^ 2 ^ 



(51) 

l (52) 
and for the momentum 



Pe 



hfd 



+ 



+ 



1 



i \d£ Z 2 + n 2 2£ 



Pv 



hfd 



+ 



n 



+ 



1 



( , 

i V dr] £ 2 + r] 2 2r\ 



(53) 



together with p^ = —ihd^. For the Hamiltonian we obtain 



2M 



-A 



E(3) 



2M 



d 2 d d 2 d 

+ 777 + TT^T + 



^ 2 + r/ 2 \9^ 2 9£ c?7y 2 t^r/ 2 d& 2 



+ 



1 a 2 



1 

2~M 



v 7 ? 2 + V 2 



(pI+p") 



+ 



ft 2 



8M£V 



(54) 
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In parabolic coordinates a potential is separable if it has the form 

u(£)+v(r]) w{4>) 



V{x) 



+ 



\/x 2 + y 2 + z 2 ' x 2 + y 2 
Prolate Spheroidal Coordinates. We consider the coordinate system 

; y/ (£ 2 — 1)(1 — vj 2 ) cos 4> = d sinh /j sin v cos (j) , /i > , < v < tt 

d sinh pi sin v sin <j) , £ > 1 , \r)\ < 1 , 

< 6 < 2tt . 



(55) 



R 



y = fV(e-l)(l 
z = -§£?? = d cosh /x cos 



r/ 2 ) sin i 



(56) 

For convenience we have also introduced the alternative representation of the coordinates in 
terms of trigonometric and hyperbolic functions via £ = cosh fj,, rj = cos v. R = 2d is the 
interfocus distance. This yields (g a b) = d 2 diag(sinh 2 fi + sin 2 v, sinh 2 /_z + sin 2 u, sinh 2 /j sin 2 u), 
and for the momentum operators we obtain 



and ; 



H I d sinh [i cosh /j 1 \ 

+■ — T2 , ■ 2 + o COth ^ ' P» 

smn /i + sin v I J 



Uld sin v cos z/ 1 

i V ai/ sinh fj, + sin z/ 2 



(57) 



-ifidfh- The Hamiltonian has the form 



"2M As(3) 

h 2 



2Md 2 
1 



1 /a 2 , d d 2 d , 

-TT2 1 T a - ? + COth A 4 -S - + "o - 9 + COt V ~ + _ 

sinh /i + sim i/ V £7> Oji ov A 



d 2 



dv ) sinh 2 /_z sin 2 z/ <9 2 </> 



2Mci 2 



1 



sinh /i + sin 2 v 



--(pI+pI) 



1 



+ 



1 



sinhV + sin 2 i/ sinh A* + sin 17 

n 2 



8Md 2 sinh /j sin z/ 



A separable potential must have the form 

u(cosh /i) + v(cos v) 



V(x) 



+ 



w(cj)) 



sinh [i + sin v sinh /_z sin i/ 



(58) 



(59) 



ci(cosh \i cos z/ + 1) gives the prolate spheroidal II coordinate system [22 



Replacing z 

|83"| , p^ , |12§| ] , which is one of the four coordinate systems which separate the Coulomb potential 
problem, c.f. the potential V^{x) in 3.2.3. Note that there is no mathematical difference between 
the prolate spheroidal and prolate spheroidal II coordinate systems. The later has a shifted origin 
such that the left focus is located at the coordinate origin. As in the elliptic coordinate system 
in two dimensions the prolate (and oblate) spheroidal coordinate system is a one parameter 
coordinate system. The prolate spheroidal coordinate has the property of separating the two- 
center Coulomb problem (Coulson and Robinson 



, Morse flTOJ and Teller JT2 
Oblate Spheroidal Coordinates. We consider the coordinate system 



x = \/(£ 2 + 1)(1 — rj 2 ) cos <j> = d cosh /j sin v cos ( 



#V(£ 2 + i)(i 



T] 2 ) sin i 



d cosh [i sin v sin < 



z = -§£77 = d sinh /j cos v 



H > , < V < 7T 

£>0 , \fj\ < 1 , 
< < 2vr 



(60) 
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(alternatively /xGlR, < i> < ir/2 pW). We again use also the alternative representation 
£ = sinh/I, 77 = cos v. R = 2d is the interfocus distance. This yields (g a b) = c?diag(cosh 2 p — 
sin 2 z/, cosh 2 p — sin 2 i>, cosh 2 /x sin 2 u), and for the momentum operators we obtain 



Pft 
and 



h d 



sinh ^ cosh /i 



1 



i I dp cosh 2 ■ : 



+ - tanh p , p p = - — + 



/i — sin" 2 
i/icl/,. The Hamiltonian has the form 



sin v cos v 1 \ , . 

+ - cot v , (61) 



i I 9f cosh 2 /i — sin 2 2 



2M 



ft 2 



2M(i 2 
1 



1 



cosh fj, — sin v \ <9/i 



d 2 , _ d d 2 d 

+ tanh fi-^ + -77377 + cot za 



2Md 2 



1 



cosh /i — sin v 



--(pI+pD 



dfj, dv 2 
1 



+ 



d 2 



dv ) cosh 2 ^ sin 2 z/ <9 2 



cosh 2 /i — sin 2 f 



+ 



1 



cosh /1 — sin v 

h 2 

8M d 2 cosh 2 p sin 2 v 

w{<t>) 



(62) 



(63) 



Here a separable potential must have the form 

u(sinh p) + v(cos v) 

= u2~^ ■ 2 - <~ u2 - „• 2 - 

cosh [i — sin z^ cosn /x sin z/ 

Paraboloidal Coordinates. The last two coordinate system are the most complicated ones and 
are similar in some of their features, c.f. (8C, 
we consider the coordinate system 



r 2 _ (m - a )(m - g)(m - a) 

a — b 

2 _ (Vi ~ b)(r} 2 - b)(r) 3 -b) 



y 



z = t;{vi + V2 + V3 - a - b) 



and |5q , 105, |106|| for further information. First 
< rji < b < 7]2 < a < r]3 

(64) 



The metric tensor is given by [P(rj) = (v ~ a )iv ~~ fy] 

'{m-m){m-m) (m - m)(m - m) (m - vi){m - va) 



(9, 



ah 



Jdiag 



The momentum operators are 

i p'(vi) \ 
4 p( m )J ' 

I P'(V2) \ 

4pm; ' 

ip' (mY 



(65) 



Pm 

Pro 



i [dm + 2 



Pm 



1(^+1 



m 




V2 


+ 


2 m 








1 




+ 


1 


1 




m 




m 


2V2 








1 




+ 


1 


1 





(66) 



Together with the dreibeins 



h mm = -1 



h mv2 = _ 1 



h V3Vs = _ 1 



P(m) 



(vi -m)(vi~V3 



P{m) 



(67) 
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we obtain for the Hamiltonian 



2M' 



2h 2 



a 



(m - m){m - m) dm 



a 



a 



dm ' {m - m){m - m) dm v dm 



a 







(m - m)(m - m) dm v dm 



(68) 



= ^ £ + ^tt)J ■ (69) 

The AVpF{m) are determined by (|To|). A potential separable in paraboloidal coordinates must 
have the form 



V{x) 



{m - m)u{m) + (m - m)v{m) + (m - m)w{m) 



{m -m){m -m){m-m) 

Ellipsoidal Coordinates. As the last coordinate system we consider 



y 



(PI 


- a 3 ) (P2 ~ a 3 ) (p 3 - 


03) 




(a 3 - ai)(a 2 - ai) 




(pi 


- a 2 )(p2 ~ a 2 )(p3 ~ 


a 2 ) 




(a 3 - a 2 )(ai - a 2 ) 




(Pi 


- ai)(p2 - ai)(P3 - 


ai) 



< ai < pi < a 2 < P2 < 03 < P3 



(ax - a 3 )(a 2 - a 3 ) 

For the metric tensor we have [-P(p) = (p — «i)(p — o 2 )(p ~~ 03)] 

/ s 1,. f(pi-P2)(pi-ps) (P2-P3)(P2-pi) (P3 -Pi)(P3 -P2) 
(Safe) = Tdiag — — , — — , — — 

4 V P(pl) P(P2) P(P3) 

The momentum operators are 



Ppl 




f d 


+ 


1 


1 




+ 


1 


1 




-tI 


\dp\ 


2 pi 




P2 


2 pi 




P3 






t d 


+ 


1 


1 




+ 


1 


1 






\dp2 


2p 2 




Pi 


2p 2 




P3 


Ppz 




f d 


+ 


1 


1 




+ 


1 


1 






\dp3 


2p 3 




Pi 


2p 3 
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Together with the dreibeins 



h PiPi = I 



P(pi) 



2 V (Pi -P2XP1 -P3 



h P2 P 2 = I 



P(P2) 



2y (P2-P1)(P2~P3) 

h p 3P s = 1 / P (P3) 

2 V (P3 - Pl)(P3 -P2) 



we obtain for the Hamiltonian 



2M' 



E(3) 



27^ 



v/npT) 



(pi - P2)(pi - P3) dpi 



P(P1) 



d 



+ 



VP{P2 







dpi (P2 - P3)(P2 - Pi) dp2 



(70) 



(71) 



(72) 



(73) 



(74) 



^(P2 



VP(P3) 







(p :1 - pj )(p :$ - p 2 ) tV-i V P ^dp 3 



dp2 
(75) 

(76) 
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The AVpp(pi) are determined by ([h]). A potential separable in ellipsoidal coordinates has the 
form of ( pTOj ) with rji replaced by pi. 

The ellipsoidal coordinate system is a two parameter coordinate system. Denoting R\ = 
02 — ai , = a 3 — °2 an d R3 = 03 — cl\ , the Ri must fulfill the constraint R\ + R%, = R% . The 
Ri describe the foci distance from the coordinate origin. By considering the limits R\2 — ► 0, 
respectively Rip 00 we can generate the seven coordinate systems described in table 1 (note 
Pi = pi - o 2 ). 

Table 1: The degenerations of the ellipsoidal coordinate system 



prolate 

spheroidal 

system 


x = (R/2)y/(£ 2 - 1)(1 - r/ 2 ) cos 4> 

y = {R/2)^{e-m-V 2 ) sin</> 
z = (R/2)t V 


i? 2 ^0 

-> i? 2 /4 


Pi -» R 2 ( V 2 - l) 
P9 — > R\ sin 2 

p 3 - RKe - 1) 


oblate 

spheroidal 

system 


x = (R/2) ^(C 2 - 1)(1 - f) 2 ) coscp 
y = (R/2) y /(e-l)(l-fj 2 ) S m ( f> 
z = {R/2)£fj 


iZ§ -> 
i? 2 -> i? 2 /4 


pi — > — i? 2 sin 2 
p 2 -(i? 2 /4)(l-fy 2 ) 


opnerical 

system 


x = r sin a cos (f> 
y = r sin 8 sin </> 
z = r cos 


/CJ — * 

#1-0 
i? 2 /i? 2 -> 


^ d2 • 2 a 
pi — * — itj sm V 

Pi —* R\ sin 2 6 sin 2 

p 3 — r 2 


Sphero-conical 
system 


x = rsn(p, fc)dn(z;, fc') 
y = rcn(/i, fc)cn(z/, A;') 
z = ?'dn(p, fc)sn(^, A:') 


R 2 ^0 
R 2 -> 
i? 2 /i? 2 - fc 2 


pi -» -fc /2 i?§cn 2 (i/, fc') 
p 2 — > fc 2 i?|cn 2 (p, fc) 
p 3 -> r 2 


Circular 
elliptic 

system 


x = (R/2) sinh p sin v 
y = (R/2) coshpcosi' 

z = z' 


-R 2 — ► 00 
R\ = i? 2 /4 


pi — » Z - it x 

p 2 — (i? 2 / 4 ) COs2 V 

P3 (# 2 /4) cosh 2 p 


Circular 

polar 

system 


x = p cos 
y = psin0 

z = z' 


i? 2 -> 00 


7\ _^ ,' 2 E>2 
pi — > Z — it! 

p 2 — > i? 2 sin 2 (j> 

P~3 ~> P 2 


Cartesian 
system 


x = x' 

y = 2/' 

z = z' 


i? 2 -> 00 
^ 00 


7\ _^ ^' 2 E>2 

pi — > z — -Kj 
p 2 — x' 2 + R\ 

,2 

P3 — y 



We must obtain the remaining three coordinate systems which are not listed in table 1. The 
paraboloidal system is obtained by the focus translation z' — > z — Ri from the ellipsoidal coor- 
dinate system by R\ 2 — > 00 together with R\/Ri — > a — b which gives three more coordinate 
systems: The parabolic system then follows from the paraboloidal system in the limit a — > b 
(the parabolic system can also be obtained from the prolate spheroidal II system). The circular 
parabolic system is obtained from the degeneration of the elliptic II system, c.f. 2.1.2. 



2.2. Separation of Coordinates in the Path Integral. 



We want to look for the coordinate systems which separate the relevant partial differential 
equations, i.e. the Hamiltonian, and, more important from our point of view, the path integral. 
In order to develop such a theory we consider according to |105| ] the Lagrangian C = 4r J2iLi hfxf 
and the Laplacian At, respectively, in the following way (where only orthogonal coordinate 
systems are taken into account) 

° 1 d (u^Mm) o \ 

A LB = 2^ 
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D 



■ £ 



a 



5i(£l, • • • ,6-1:6+1: • • ■ :£d)/(6 



96 



(77) 



necessary for the separation []7|, 105]. We introduce the Stackel-determinant (7J, 101, [L0 



where {£} denotes the set of variables (£i, • • • ,6d), an d the existence of the functions /», <fc is 



(78) 



D 



S{{i}) = det(^) = [] 



htm) 



i—l /i(6) 
and abbreviate Tj = f[/fi- Then 



i+1, 



55 £({£}) 



i(Ci> • • • ,6-i)6+i: • • ■ >£d) = Mi(Ci> • • • ,6-1)6+1: • • ■ >£d) n 



0) ' 



(79) 



3=1 

«^3 



which fixes the functions gfj. Introducing the (new) momentum operators = y(c?6 + hfl/fi) 
we obtain according to the general theory the following identity in the path integral |5^] by 
means of the space-time transformation technique (Duru |28j| , Fischer et al. |^0[, Refs. 
Kleinert 85], and Pak and Sokmen [ 110 1 ) 

&(*")=«" 



D 

n 



hiT>£i(t)exp 



D &(*")=«' 



M 



E^I-AVpi,({e}) 



i=i 



dt 



n / ^B6w«p u' 



(*')=^ 



eft 



*— 1 Ctl\\— f 



[ i M D f 2 h 2 u 

K L 1=1 i=l 



2 D 



ds 



(80) 



Here we have adopted the following lattice formulation of the path integral (c.f. DeWitt [ pBj| , 
Feynman |37[|, Gervais and Jevicki [47[|, Refs. (5^, |6l], McLaughlin and Schulman |J92|| , Mizrahi 



and Omote [107] for alternative lattice definitions) 

q{t")=q" 



K(q",q';T) 



q{t')=q 



J \/g V PF q (t) exp | ^ ^ 

/ M \ND/2N_l 



Y h ac(q)h cb (q)q a q b ~ V(q) - AV PF (q] 



3=1 L 



3=1 

M. 



—hU&Kdqj-JAqfAq? - eV(qj) - eAV PF ($j) 
Ze 



(81) 



This path integral formulation will be used in the sequel. 
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3. Path Integral Formulation of Smorodinsky-Winternitz Potentials. 
3.1. Two-Dimensional Smorodinsky-Winternitz Potentials. 

In this subsection we discuss the four two-dimensional Smorodinsky-Winternitz potentials 
They are characterized by having three functionally independent integrals of motion, i.e. there 
are two more operators related to these integrals of motion which commute with the Hamiltonian. 

We will state if ever possible both the propagator and its spectral expansion into the wave- 
functions and energy-levels, x denotes a two-dimensional variable: x = (x,y) = (xi,X2), and 
without loss of generality it will be assumed that x 6 IR 2 ,IR 3 . 

In the table 2 we list the two-dimensional maximally super-integrable Smorodinsky-Winternitz 
potentials together with the separating coordinate systems f42fl . The cases where an explicit path 
integration is possible are underlined . 

Table 2: The two-dimensional maximally super-integrable potentials 



Potential V(x, y) 


Coordinate 
System 


^V + /) + 2 * M (V + V) 


Cartesian 

Polar 
Elliptic 




Cartesian 


Parabolic 


y « , 1 ( k l-\ , k i-\ ) 

y/x 2 +y 2 4M ^ x 2 + y 2 \ \/x 2 + y 2 + x \/x 2 + y 2 — x J 


Polar 

Elliptic II 
Parabolic 


a [2/ 4> (j>\ 
^4 = -- + y-( v /3icos- + /3 2 sin-j 


Mutually 
Parabolic 



3.1.1. We consider the potential (kip > 0) 

M 



(82) 



This potential is separable in three coordinate systems. We obtain the path integral formulations, 
c.f. ( |8l|) for their lattice definitions 

K^\x",x';T) 
Cartesian Coordinates: 



x(t")=x" 



T>x(t) exp 



x(t')=x' 



M ,-_2 
~2 



2 1 1.2 _ 1 



2-2* ft 4 1 ^2 4 
(X — U! X I — — — q H s-^ 

V ' 2M\ x\ x\ 



dt 



(83) 
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Polar Coordinates: 

P (t")=p" 4>{t")=<i>" 

pVp(t) 



v<t>{t) 



p(t')=p' 



xexp - 



Ht')=<t>' 



M 
~2 



I -2 , 2 12 2 2\ 
— [ft + p <j) -Up) 



h 2 



2Mp 2 \ cos' 



1 h.2 1 

4 i ^2 4 

i "r . o , 



sin 



4 



(84) 



Elliptic Coordinates: 



?(*")=?" 



r?(*"W 

/ 

r,(t')=v' 
M 



£(*')=«' 



xexp - 



y X>£(i) y Pr/(t)d 2 (sinh 2 ^ + sin 2 r ? ) 



— d 2 ^(sinh 2 £ + sin 2 i])^ 2 + rj 2 ) — w 2 (cosh 2 £ cos 2 r] + sinh 2 £ sin 2 rj 



ti 2 



2 1 
4 



+ 



i-2 _ 1 
ft 2 4 



IR 



dE 



2Md? \ cosh 2 £ cos 2 r\ sinh 2 £ sin 2 77 



eft 



(85) 



-iET/h 



ds" 



€(o)= 



xexp - 



2M 



— ((£ 2 + ?) 2 ) — w 2 (cosh 2 £ sinh 2 £ + sin 2 77 cos 2 rj) ) + i^cf 2 (sinh 2 £ + sin 2 77) 



n(o)=v' 

,2/„„„u2 



Vr)(s) 



(k 2 



2 i> 

4' 



1 



1 



cos"' 77 cosh £ 



+ (fc 2 



+ 



sin 2 77 sinh 2 £ 



(is 



(86) 



Note that in the last line a pure time-transformation [84] has been performed. 

We are going to discuss the path integral (Q) in some detail in order to demonstrate the 
relevant technique. All subsequent path integral evaluations are similarly done and our method 
remains the same whether polar, parabolic, or other coordinate systems are used. 

The path integral solution in cartesian coordinates is straightforward by applying the path 
integral solution for the radial harmonic oscillator for the coordinates X12, respectively, c.f. ( p^ , 
|93[ ). In the case of polar coordinates one first separates the 0-path integration by means of the 
path integral identity of the Poschl- Teller potential (0 < < vr/2) || 28, 72] 



<p{t")=r 

I 

<t>(t')=<t>' 



V<f>(t) exp I - 



M 



h 2 ( a 



S 1 



2M V sin 2 , 



cos^ 



dt 



(87) 



ra=0 



and $>^'^ ((/>) denote the normalized Poschl- Teller wave-functions 



2{a + (3 + 21 + 1) 



UT(a + (3 + 1 + 1) 



1/2 



T(a + l + l)T((3 + l + l) 
x (sin 0) Q+1/2 (cos <j,)P+y*pM (cos 2(f)) 



The Pn (x) are Jacobi polynomials ( fl49fl , p. 1035). Note that the genuine Poschl- Teller poten- 
tials requires a, (3 > \. Otherwise one can call it an attractive Poschl- Teller-like potential. This 
rule of separation of variables in the path integral consists of the following path integral identity 
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© (g = T\g 2 , {9ih = 9 = m) 



x(t")=x" 



r(z)^g(z)Vz(t) J Vx(t) 

z(t')=z' x(t')=x' 

' ft" M 



x exp 



[f(z')f(z")]- d/2 J dE x *W)* x ( 



V{x, 



x exp 



z(t")=z" 
z(t')=z" 

Ex 



W(z) + AW(z) 



dt 



^-(g.-z) 2 -W(z)-AW(z) f2{: . 



dt 



where the x-path integration is given by 

x(t")=x" 



x(t')=x' 



Vx{t) exp < — 

n Jt 



yi 2 - V(x' 



dt 



dE x <£> x (x')$ x (x")e 



(90) 



The remaining /5-path integration ( |97| ) is again of the radial harmonic oscillator type. The radial 
l/p 2 -dependence must be incorporated into the path integral in terms of a functional weight 
fJ>x[p 2 ] ( |98| ) in order to guarantee a proper definition of the corresponding lattice formulation of 
the radial path integral |119| . As it turns out, this Besselian path integral formulation jTfj] 
with the functional weight fJ-x[p 2 ] f° r radial path integrals can be interpreted as an additional 
potential term, and vice versa J|(| |8(| . We also want to point out that a similar reasoning must 



be done in the case of the Poschl- Teller and modified Poschl- Teller potential, respectively [41 



Hence we apply the radial harmonic oscillator the path integral solution and its expansion into 
wave- functions according to (A an arbitrary parameter |2^, 48, 113 , |119| ) 



r(t")=r' 



T>r(t) exp 



r{t')=r' 



1 



M 
Y 



—f 2 - h 2 



A 



M 



±-^r 2 )dt 
2Mr 2 2 / 



r(t")=r" 

J nx[r 2 ]T>r(t)eyip 

r(t')=r' 

Mojy/r'r" 



\M 
~2h 



(r 2 -iv 2 r 2 \dt 



ih sin ioT 
2Mlo 



exp 



h 

x exp 



V r'r" e 



- — r /2 +r" 2 cot 
2in 

AwT(2n+X+l) n - 



h 



Mujr'r" 
i?isin loT 



(91) 
(92) 



n=0 

Mlo 

~ ~2h 



n\ {Muj . 

r 7 

r(n + A + l)V h 

(r' 2 +/))4 A) (x^i A) ( 



A 



Muj //2 



(93) 



The expansion into the wave-functions has been obtained by means of the Hille-Hardy formula 
(H, p.1038) 



r"/ 2 
l - 1 



exp 



1 1 + t 

2 V yJ l-t 



I, 



°° fn n \ e -^(x+y) 

E r( - I ■ ^ xyTl 2 L n {x)L n {y) . (94) 



2^/xyi 

l-t J jL ^ n T(n + ot+lY 



The L n (x) are Laguerre polynomials p. 1037). We therefore obtain 
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K^\x",x';T) 

Cartesian Coordinates [29, E8|, 113, 119]: 



Moj 



\Tt sin ojT 



2 2 



oo 

= E e 

m,n2=0 
Polar Coordinates: 



i=i 

iE N T/h^ , ../ .,/ 



+ xf ) cot uT 



\i?isinu>T 



ni,ri2 { x li x 2)^ 'ni,ri2 i x l j ^2) 



p(i")=p" 



1 



_^ $ tfcfc 3 ,±fc 1 )^/ ) ^±fe,,±fc 1 ) (0 // ) / p/0(t)exp 



n=0 



p{V)=pf 



M 



\ 2 1 
2 A ~ 4 



2 (^- w V)-fi 



p(t")=p" 



(p'/)- 1/2 E^ ±fc2 ' ±fcl) (^)4 ±fe 



n=0 

00 



/_i A [p 2 ]Pp(t)exp 



P (t')=p' 



iM /■* 



2/i 



{p 2 -uj 2 p 2 )dt 



ihsinuT 



J2 $tfck a ,±fc 1 )^/)$ttfc a ,±ki)(0//) exp 



n=0 



— (p +p Jcot^T 



/ Mup'p" 
\ih sin wT 



m,n=0 



V) , 



(95) 
(96) 

dij 

(97) 
(98) 

(99) 
(100) 



where A = 2n ± fc^ ± &2 + 1. The wave- functions and the energy-spectra are given by 



Cartesian Coordinates: 

2Mu A (Muj\ ±ki ' 2 



i=l 



E N = Tluj{2N ±k 1 ±k 2 + 2) 
Polar Coordinates: 

Wp,0) = 4 ±fc2 ' ±fcl) (<^ 



E N = hu(2N ±k 1 ±k 2 + 2) 



y T(ni ±ki + l' 
N = m + n 2 , 



l/2±fcj 
x i exp 



2h V n ' { h 



(101) 
(102) 



2Muj/Muj 2 \ A / 2 

—It p 



m! 



T(m + A + 1) 



exp 



Muj o 
2T" 



4J>(^fV 



N = m + n , 



(103) 
(104) 



and <J?n a '^ (</>) of (^) with a = ifo, /3 = i/ci. Let us remark that the wave-functions have been 
normalized in the domain [0, 00). The positive sign at the k{ has to be taken whenever k, L > ^, 
i.e. the radial potential term is repulsive at the origin, and the motion takes only place in the 
domain X{ > 0. If < \k%\ < 5, i.e. the radial potential term is attractive at the origin, both 
the positive and the negative sign must be taken into account in the solution. This is indicated 
by the notion rbfcj in the formulas. It has also the consequence that for each ki the motion can 
take place on the entire real line. In the present case this means that in cartesian coordinates 
we must distinguish four cases: i) xi,x 2 > 0, ii) x% > 0, x 2 £ 1R, hi) x\ 6 IR, x 2 > and 
iv) (x\,x 2 ) E IR 2 . In polar coordinates the same feature is recovered by the observation that 
the Poschl- Teller barriers are absent for < |. We will keep this notion in the sequel for all 
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following Smorodinsky-Winternitz potentials. 
3.1.2. We consider the potential {k\^ > 0) 



V 2 (x) = f u 2 (Ax 2 + y 2 ) + M + 



(105) 



which is separable in two coordinate systems. For k\ = this potential is also called Holt- 
potential [66]. We have 

K^\x",x';T) 
Cartesian Coordinates: 



x(t")=3" 

Vx{t) exp < - 

x{t')=x' ^ 

Parabolic Coordinates: 



M rJ> 



x — u> 2 (4:X 2 + y 2 ) ) — k\x 



h 2 k\ I 

2M y 2 



dt 



(106) 



v(t")=v" 

I 

v(t')=v' 



€(*")=£" 



| V V (t) J Vi{t){e + rf 



expU/, f ((e 2 + ^ 2 )(e 2 + ^ 2 )-- 2 ((e 2 



1 /v2 



k 2 - I 
2 ^2 4 

2Mr) 2 i 2 



(107) 



JR 



-iET/h / dg 



r?(0)=77' 



«(0)=?' 



xexp - 



^((e 2 + ? ? 2 )-^(e 6 + r ? 6 ) 



ft? ( ^2 4 , ^2 I 



(is 



(108) 



In the x-path integration one makes use of the shifted harmonic oscillator path integral solution 
[37, |38| , and in the y-path integration one makes use of the radial harmonic oscillator path 
integral. The expansion into the wave- functions in x is obtained by means of the Mehler formula 
(Hi, P-272) 



-(x 2 +y 2 )/2 



E 

n=0 



1 (z 



re! \ 2 



H n (x)H n (y) 



: exp 



Axyz - (x 2 + y 2 )(l + z 1 



2(1 -z 2 ) 



(109) 



The H n (x) are Hermite polynomials (||49j, p. 1033). We have for the harmonic oscillator path 
integral the kernel and its expansion into wave-functions according to (we omit the Maslow 
indices) 



x(t")=x" 



Vx{t) exp 



x(t')=x' 



2Tt\h sina;T 



exp 



Moj 



[x' 2 + x" 2 ) cot ujT — 2 — 

sinw 



n=0 



-iwT(n+l/2) / Mw xl/2 



2 n re! V 7T^ / 



a/ I 



x exp 



(110) 
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The £ 4 , 7j 4 , £ 6 , ry 6 -terms make the path integral in parabolic coordinates intractable. For k\ = 
the anharmonic sextic potential problems in £ and rj can only be treated for a zero separation 
constant, c.f. the appendix. The result in cartesian coordinates has the form (x = x + fci/4Mw 2 ) 
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K^\x",x';T) 



Mu \ x l 2 



2ir\h sin toT 



Muj 



2\Tism.ujT 



exp 



exp 



\Mu) 



i~l2 i ~ll2\ rr r,~l~ir 

[X +X ) COS UJl — IX X 



2ft, sin wT 



Mid, ,2 „2, „ 

— (y +y JcoscjT 



fc2 V iti sin u;T 



ni ,ri2=0 



(112) 
(113) 



with the wave-functions and the energy-spectrum given by 



1 2Muj (Muj\ ±k2 / 2 



n 2 \ 



I2Muj 1 



T(n 2 ± k 2 + 1 
eXP l" — 



Muj r 

~2~h ] 



r(±fc 2 ) 

ri2 



£7, 



711,712 



^ V ^ 2 n ini! 
fiw(m + 2n 2 ± fc 2 + I) + " 1 



x 2 \H, 



I2Muj 



n 



k 2 

2) + 8Muj 2 ' 



(114) 
(115) 



Note the cases whether k 2 > \, respectively < k 2 < \. 
3.1.3. We consider the potential (ki i2 > 0) 



V 3 (x) 



a 



+ 



ti 2 



p _ l 



+ 



k 2 - 1 
K 2 4 



\/ x 2 + y 2 4M ^Jx 2 + y 2 \ \/ x 2 + y 2 + x \/ x 2 + y 2 — x 



(116) 



This potential is separable in polar, parabolic and elliptic II coordinates. A path integral 
discussion for the pure Coulomb case is due to Inomata [68] and Duru and Kleinert [^TJ. We 
have 

K^\x\x'-T) 

Polar Coordinates: 

P {t")= P " 



pDp{t) 



P (t>)=p> 

x exp 



<f>(t')=4>' 



h 2 k 2 



M . 2 9- 2 . a 
— (p + p </>) + - 
i p am p \ cos ^ r gin 



2 _ 1 1.2 1 n 
4 + ^2 I _ 



(117) 



Elliptic II Coordinates [95]: 

e(*")=e" 7?(t")=77" 

= /" Pf(t) J P7 ? (t)d 2 (sinh 2 e + sin 2 ? ? ) 

77(t')=77' 

^(sinh 2 £ 2 + sin 2 ?y)(£ 2 + ?) 2 ) + 



x exp 



ft Jt' 



a 



ft 2 



d(cosh £ + cos rj) 

2> 



' (k\ + k\ - \) + {k\ - k\) cos j] 



4Md(cosh 2 £ - cos 2 rf) 
+ 

Parabolic Coordinates [Bill: 



I sin 2 7] 

(k 2 + k 2 -\) + {k 2 -k 2 ) cosh ^ 



sinh 2 £ 
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ri(t")=r)" 



v(t')=v' 
x exp 



e(*")=r' 

/ T>r,(t) J 2>£(t)(£ 2 + r/ 2 

e(f)=c 



M ,,2 2W2 -2n 2 " & /^l- 



1 1,2 1 ' 
1 , K 2 4 



+ 



(2ij . 
(119) 



If the domain of the angular variable is ^ 6 (0, 2tt) there are two separated domains (half- 
planes) in the case of k\ > \. A rotation of the axis about it/ 2 interchanges the roles of 
k\ and ki- In order to evaluate the two path integrals in polar and parabolic coordinates a 
space-time transformation must be performed. The result in polar coordinates has the form 
[k = a\j 1 —M /2E /h, A = m+ (1 rb k% ± &2)/2, only the radial Green function can be explicitly 
evaluated] 



G^\x",x';E) = i / dTe iTE / n K^\x\x'-T) 
n Jo 

= i( r y')-V2 J $ ttfc2,±fe0 $ (±* 2 ,±fei) ^ 



z=o 



p(t")=p" 



dT jTE/K 



Vp(t) exp 



P(t')=p' 



i /•*" / M 

h 



2 P p 2Afp2 



oo I oo 



$ (±fc 2 ,±fc x ) $ (±fc 2 ,±fei) 



m r(| + A - k) 
2~E T(2A + 1) 



W K , X ( v 7 ^^^ ^ ) AT- A ( ^8ME ^ 



y . y ^n, m (^ / ,^ / )K,m(p // ,0 // ) + 



n=0 I m=0 



£",71, — -E 



ft 

r P h 2 p 2 /2M-E 



(120) 
(121) 



In the c/>-path integration the path integral solution of the Poschl- Teller potential has been used. 
In the radial path integration the path integral solution for the radial Coulomb potential was 
used (Chetouani and Hammann |l{|, Duru and Kleinert (3^, 31], and Refs. |56|, |11S| 1), and the 
Green function is expanded according to (Kratzer potential, k = ct\J — M/2E /%) 



x(t")=x" 
r 

dT jTE/K / Vx{t) 

exp 

x(t')=x> 

i / m r(| + x — k 

oo 



M . 2 a 

~2 X + Jx\ ~ 2M~c 2 ~ 



h 2 A 2 -i 



\dt 



2E T(2A + 1) 



■w K 



h 



V-8ME ) M K \ ( V-8ME ^ 



Ti 



n=0 



E n — E 



JR 



E p — E 



(122) 
(123) 



x <5> denotes the smaller /larger of x',x". The W K; \(x) and M Kj ^(a;) are Whittaker functions 
d|49t| , p. 1059). The bound state wave- functions and the energy-spectrum of the Kratzer potential 
are given by (a = h 2 /Ma) 



1 



n + A + i 



nl/2 



aT(n + 2A + 1) 
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2x 



En 



a(n + A + \) 
Ma 2 



A+l/2 



cxp 



X 



a(n + 1) 



2x 



a(n + 1) 



and for the continuous spectrum with E p = p 2 h 2 /2M one has 



yjx) = v 2 , —^-ex-p — )M i/avX (-2ipx) . 

py 1 v&(2A + l) \2apJ 1/ap ' M F ' 



(124) 
(125) 

(126) 



Therefore the polar coordinate wave-functions and the energy-spectrum of our potential problem 
are given by (a = M/ah 2 ): 



2-l/2$(±*a,±ki) ( 



\2 



2m! 



2p 



a 2 (m + A + i) 3 r(m + 2A + 1) 
A 



1/2 



a(m + A + 5) 
Ma 2 



cxp 



a(m + A + |) ' 



L (2A) 



2p 



a(m + A + I) 



2fi, 2 (m + A + 



l2^J Mi ^A(-2w) • 



, (127) 
(128) 
(129) 



The $n ±fc2 ' =l:fcl ^((/ , /2) are the Poschl- Teller wave-functions (|3£|). 

In parabolic coordinates we obtain by performing a time transformation and applying the 
path integral solution of the radial harmonic oscillator (92) in the £ and 77 variable, respectively 

(w = y/-2E/M) 



i TOO 

^ / dTj TE / h K^(x",x';T) 
n Jo 



V (s")=v" 



{(«")=*" 



Dr/(s) 



2?£(s) 



x exp 
2 



r?(0)=7?' 



M 



C(o)=e 



fr 2 ( k\ 4 ^2 I 



+ 



sin 2 ws" 



Jias" /h 



x exp 
00 ^ 



E 

ni,n 2 =0 



n.1,712 



/ 2 e //2 ,2 ,/2s , ,/ 

(t',T/)* nitna {Z",r/') , /-oo 
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h 2 



Mujrj'r] 1 ' 
\h sin ws" 



K — E 

- t ^ni,ri2 ■*- / 



+ dp dC— — - — , 
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(130) 



where are (N = n\ + n. 2 + ± ^2) + 1, a = h 2 /Ma) 



ni!ra 2 ! 



a 2 N 3 r(m ± fci + i)r(n 2 ± fe 2 + 1) 



1/2 



E 



x 

Ma 2 
h 2 N 2 



(-L) ^ (JL) exp ( - ?±£\LW ( * 
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X M -i(l/a+C)/2p,±fci/2(- i K 2 )^-i(l/a-C)/2p,±fc2/2(- i W 2 ) , 

^ ft 2 2 

^=2M^ • 



(133) 
(134) 



In order to extract the discrete spectrum one applies (c.f. [20, pq| ) the Mehler formula [E 
For the continuous spectrum one applies the following dispersion relation (c.f. ||, p. 158, £ 
p.414, H, p.884) 

1 



ana 



exp 



(x + y) cot a 



2^ 



ana 



1 r r(| + A* + ip)r(i + M -ip) 



2n^xyM r 2 (l + 2//) 



e -2ap+^ M+ . p iM (_2ix)M_ ip , M (+2iy)dp . (135) 



and performs the variable substitution (^1,^2) ► {PiiPn) *— > [(V a + C)/2p> (1/a — C)/2p]- C is 
the parabolic separation constant. In elliptic II coordinates ET^ 3 ) is not explicitly soluble and 
we obtain only the path integral identity 



dE 



C(«")=f" 



«7( S ")=r?" 



■n(o)= v > 



xexp - 



— (£ 2 + rj 2 ) + Ed 2 (cosh 2 £ — cos 2 rj) + ad(cosh £ + cos rf) 



K 1 (kj -{)(!- cos rj) + (fc 2 2 -{)(! + cosr?) 
2M I 



sin 2 77 



"l ~ 4 



+ - 

No further evaluation is possible. 

3.1.4- We consider the potential ({h,2 £ It) 



)(1 - cosh£) + (A; 2 - i)(l + coshf)' 



sinh £ 
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+ \ / - ( Pi cos — + P 2 sin • 
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This potential is only separable in mutually orthogonal parabolic coordinates and again a time 
transformation must be performed. We obtain [to = \J —1EjM , only the Green function can be 
explicitly evaluated, (£, fj) = (£ — 0i/E, rj — 02 /E)] 
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dT jTE/K 



Vx(t) 



xexp - 
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x exp 



f(e 2 + , 2 )(e 2 + , 2 ) + 2 a - ( ^ + /327?) 



dt 
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xexp - 
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(139) 
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h 2 p 2 /2M -E 

and S e o denotes the summation over even and odd states, respectively; the bound state energy- 
levels are determined by (N = n\ + n 2 + 1, n\, n 2 G INq, ojn = \/—2En/M) 



3 2a , n /? 2 + /? 2 2 



. 



(142) 



For a > the discriminante of this cubic equation shows that there is one real root ujn for each 
N such that the bound state energy-levels En are given by 
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The corresponding bound-state wave-functions have the form 
I'm ,n 2 (£,</■) = 
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iV = n\ + ?i2 must be chosen in such a way that N = n\ + n 2 is an even number [130]. The 

n(0) ,'„-i7r/4 



continuous functions ^i 6 /-*^ (£, ??) are given by [a = H 2 /M(a - m(f3 2 + /3|)/^V)] 
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For the determination of the continuous spectrum one uses the dispersion relation ([|49|, p. 896) 
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and performs the variable substitution (^1,^2) l—> (p£,Pr)) l— > [(V^ + C)/2p> — C)/2p]- C is 
the parabolic separation constant. The D u (z) are parabolic cylinder functions, and the E^\z) 
and E^\z) are even and odd parabolic cylinder functions ||, respectively. 

Note that in the evaluation of the path integral one has to take into account that by using the 
linearly shifted harmonic oscillator solution for the £- and 77-variable, respectively, one actually 
uses a double covering of the original (x,y) 6 M 2 -plane, i.e. u = (£,fj) € IR 2 . Furthermore we 
have taken into account that our mapping is of the "square-root" type which gives rise to a sign 
ambiguity. "Thus, if one considers all paths in the complex z = x + iy-plane from z' to z" , they 
will be mapped into two different classes of paths in the u-plane: Those which go from u' to u" 
and those going from v! to — u". In the cut complex z-plane for the function \u\ = \/\z\ these 
are the paths passing an even or odd number of times through the square root from \z\ = and 
\z\ = —00. We may choose the u' corresponding to the initial z' to lie on the first sheet i.e. in 
the right half u-plane). The final u" can be in the right as well as the left half-plane and all 
paths on the z-plane go over into paths from u' to u" and those from u' to — u" " |3l|| . Thus the 
two contributions arise in (140). 



3.2. Three-Dimensional Maximally Super-Integrable Smorodinsky-Winternitz Po- 
tentials. 

We discuss in this subsection the five three-dimensional maximally super-integrable potentials. 
They are characterized by having five functionally independent integrals of motion. In the sequel 
x denotes a three-dimensional coordinate: x = (x,y,z) = (xi,X2,xs). 

In table 3 we list the three-dimensional maximally super-integrable Smorodinsky-Winternitz 
potentials together with the separating coordinate systems (where the italized coordinates sys- 
tems were not mentioned in p5| ). The cases where an explicit path integration is possible are 
underlined. 



3.2.1. We consider the three-dimensional potential (^1,2,3 > 0) 



Vx{x) = -ux +— - 
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+ 



+ 
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Note furthermore that the pure three-dimensional harmonic oscillator [73] has the same sepa- 
rating coordinate systems as V\{x) and therefore the additional centrifugal terms do not spoil 
this property. 

We formulate the corresponding path integral in the coordinate systems which separate the 
path integral for this potential problem. We have 

K^\x",x»-T) 
Cartesian Coordinates: 
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Table 3: The three-dimensional maximally super-integrable potentials 



Potential V(x, y, z) 



Coordinate System 



M 



Vi = '~-to 2 x 2 + -— 
2 2M 



t,2 / h.2 1 pi h 2 — 1 



+ 



+ 



Cartesian 

Spherical 
Circular Polar 
Circular Elliptic 
Conical 

Oblate Spheroidal 
Prolate Spheroidal 
Ellipsoidal 



w = 0, fe| — 1/4 = 
u = 0, k\ -1/4 = 



Paraboloidal 
Parabolic 
Circular Parabolic 



M 



fi 2 ( k 2 \ k\ | 



+ 



x^ 
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Circular Polar 
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+ 
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'h.2 1 h 2 1' 
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x^ 
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+ 



/Co 



+ 
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h 2 ( 
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2M v / x 2 +y 2 y Vx 2 +7 + a; V^ 2 + V 2 
01 = 5(^2 + ^1 + 3); /5| = 5(^2 — ^1 + 3) 
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Circular Elliptic II 

Circular Parabolic 
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Circular Polar Coordinates: 
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Circular Elliptic Coordinates: 
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2M<i 2 y cosh 2 sin 2 v V cos 2 </> sin 2 <j> 4/ sinh 2 cos 2 v / 

oo 

(d 2 cosh//' cosh /2" sin sin P"y x l 2 ^ k2 ' ±kl) {(j)')^ k2 ^ kl) {f) 



dt 



In 



Ar 2vr/i 



ra=0 

p,{s")=p," P(s")=P" 



poo r r 

J ds" J Vfl(s) J VP(s) 



P(0)=P' 



(152) 



(153) 



(154) 



f i f lifi 

exp < — J -^{{p 2 + v 2 ) — w 2 ( cosh 2 /i sinh 2 p, + sin 2 Pcos 2 P)^j + E'^cosh 2 p, — sin 2 P) 



2M 



(A 2 



2 1- 



1 



1 



cos 2 sinh /i 



Prolate Spheroidal Coordinates: 



**(*")=**" 

/ 

Mf)=A»' 



x exp 



i/(t")=i/" 

/ 

u(t')=v' 

M 



sin 2 z/ cosh 2 fi 



Ht")=4>" 



ds \ (155) 



J T>n(t) j Vv (t)d 3 (sinh 2 n + sin 2 i/) sinh fj, sin i/ ^ P^(t) 



*(f)= 



t' 



<i 2 ( (sinh 2 fi + sin 2 z/) (/i 2 + P 2 ) + sinh 2 // sin 2 z/(/> 2 
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0/ i0 0\ 

-co (sinh ^, sin Z/ + cosh p cos 



Ti 2 



1 



/P - i 



1.2 1 
4 , K 2 



+ 



2Md 2 I sinh 2 // sin 2 v \ cos 2 sin 2 4 
Ellipsoidal Coordinates: 

pi(t")=Pi' p 2 (t")=P' 2 ' P3(*")=P 3 ' 



+ 



t 2 1 

ft 3 4 



Ppi(t) 



Vpz{t) 



cosh 2 /Lt cos 2 f 

(P2-Pl)(P3-p2)(p3-Pl) 



Pi(i')=Pi 



P2(t')=p' 2 

3 „ f 



P3(*')=P3 



M . 2 l_ 



8 v /-P(p 1 )P( / 9 2 )P(p 3 ) 
/Mw 2 n/ , ?i 2 



(it 



(156) 



dtj . 

(157) 



In this example we have explicitly written down the separated (£, ry)-path integrations in the 
case of the oblate spheroidal coordinates. The case of the prolate spheroidal coordinates is, of 
course, similar. For the separation in the ellipsoidal coordinates A{p) denotes (ojfc = o« — a^j 



k 2 - - 

Mp) = «3i«2i + 012^32 

p-ai p- a 2 



k 2 - I 

ft 2 4 i ^3 4 



+ 013^23" 



/l " ■< 



p-a 3 



(158) 



The separation of variables then is performed by means of (pO|). In the case of the conical 
coordinates the separation of variables in the path integral becomes obvious if we consider for 
fixed R the identity 



a(t")=a" f}(t")=(3" 
ot{t')=a' P{t')=P' 



Vf3(t){k 2 cn 2 a + k' 2 cn 2 p) 



xexp 



-^-R 2 (k 2 cn 2 a + k> 2 cn 2 p)(<r + f) 



K 1 



h 2 1 h 2 1 

ft l 4 _|_ ft 2 4 



+ 



k 2 - i 

ft 3 4 



IR 



dE 
2nh i 



2MB 2 ^sn^dn 2 ^ cn 2 acn 2 /3 dn 2 asn 2 /3 / 
a(s")=a" P(s")=P" 



dt 



-iET/h 



ds ' 



a(0)=a' 



x exp < - 

n- Jo 



Va(s) 



h 1 



Vf3{s) 



1 



+(fc 2 



/ A; 2 k 
' + 



i2 



\cn 2 /3 cn 2 a 



sn 2 a dn (3 
1 



k 



sn 2 /3 dn' 



a 



ds 



(159) 



We find that the path integral for this potential can be explicitly evaluated in the cartesian, 
circular polar and spherical coordinate system [jlC|J. The others cannot be explicitly evaluated, 
but are however connected through the above path integral identities. We obtain (c.f. for path 
integral discussions for special cases of ^1,2,3 e.g. Carpio-Bernido and Bernido [[l(], |i~2| ]) 

K^\x",x'-T) 
Cartesian Coordinates: 
( Muj ^ 3 3 



ih sin loT 



II V x 'i x 'l ex P 



4 = 1 



^(xf + x- /2 )cot ujT 



( Mujx^x'I 
Vi^sinwT 



(160) 
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n\,ri2, n 3 =0 

Circular Polar Coordinates {X\ = 2n ± &i ± k2 + 1): 



(161) 



i/isinwT 



cxp 



2ih 



z + z ) cot wT 



Muz'z" 
ih sin cjT 



x J2 $ { ± k2 ' ±kl) (^')$i ±fc2 ' ±fcl) (</>") exp 

n=0 



— [p +p ) cot uT 



f Mujp'p" 
^ \ih sinwT 



= £ e-^ r / R * n , tnjn ,(^y,y)* n , m , n ,(^ / ,/,o 

n,m,n z =0 

Spherical Coordinates (A2 = 2m + Ai ± £3 + 1): 

00 00 

= (rV / sm0 , S me ,/ )- 1/2 E $ n ±fc2 ' ±fcl) (0 /, )^ ±fc2 ' ±A:i) (0 / ) £ ^'^(fl'')*^'** 8 ^) 

n=0 m=0 

MrV" 



(162) 
(163) 



Mu) 

x — — exp 

msm ujT 



2ih 



(/ +r" )cotuT 



I\ 2 



ih sin loT 



(164) 
(165) 



l,n,m=0 



with the wave-functions and the energy-spectra given by (N is the principal quantum number) 
Cartesian Coordinates: 

/2MuA 3 / 2 A /Moj\ ±k >/ 2 



9 



ni,n 2 ,n 3 



(xi,X 2 ,X 3 ) 



1—; n(- 



r(raj ± fej + 1' 



2h 



\ n 



i=l 



f?jv = ?M2iV + £(1 ±fc ») ) ' A f = ni + n 2 + n 3 , 



(166) 
(167) 



Circular Polar Coordinates (Ai = 2n ± fci ± hi + 1): 

™ n 



,m,n z \ 



n- 



T(n z ±k 3 + l 



_ z l/2±k 3 exp 
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I2Muj (Moj^\ x ^ 2 



T m + Ai + 1 



V h V h 

E N = huj(2N + J2(l ± fei) ) , iV = m + ra + 



exp 



I2T 



i=l 



Spherical Coordinates (A2 = 2m + Ai ± £3 + 1): 



2Mw/Mw 2 \ X2 ' 2 

\ h \ h 

E N = ftw(W + £(l±^ 



■ exp 



Muj 

~2~h 1 



V r(z + A 2 + i) 

, N = n + m + I , 

and <J?i ±fc2 ' =l=fel ' l ((/>) and Qm 1 '^ (0) denote the wave-functions (|88|), c.f. Calogero |]] and Evans 



i=i 



(169) 



L^(^r 2 ), (170) 
(171) 
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1 35, 36] for the Schrodinger approach. Note that our solution in circular polar coordinates is 
new. The L>-dimensional generalization of this potential is also maximally super-integrablc 
and it is immediately obvious that it is at least separable (and exactly soluble) in cartesian and 
spherical coordinates (and any combination of circular spherical subsystems). 

Let us furthermore consider the two following special cases of the potential V\{x) where 
ijj = 0. Here we have a continuous spectrum instead of a discrete one. 

(i) k 2 - \ = 0: We have (A = 2n ± ki ± k 2 + 1) 

K^\x",x'-T) 
Paraboloidal Coordinates: 



vs(t")=v% 



x cxp 



V m (t) 

• 3 

is 



M 



v m (t) J v m (t) 

V3(t')=v' 3 

ft 2 B(th) 



8^/-P( m )P( m )P( m ) 



Parabolic Coordinates: 

V (t")=v" e(t")=e" 
I T>n(t) J vat)(f + v 2 )^ I v<t>{t) 



AVpF(r]i 



dt 



(172) 



fi 2 / fcjf 



2M£ 2 rf \ cos- 



2_I p _ 1 



sin 



(173) 



(^'wr l/2 E $ i ±fe - ±y (^e 



2,±fcl)/ ///X 



n=0 



JR. 



x exp 



2vrn ( 
ft 



r?(s"W 



i{s")=i" 



-LET/h / ds « 



£>r?(s) 



V(0)=rj> 



e(o)=€ # 



f f((e 2 +f)+^ 2 +r/ 2 )-ft 



2 1 



1 1 

2M V^ + ^ 



(is 



n=0 



IR JO 



(174) 
(175) 



In paraboloidal coordinates AVpjr^j) are determined by ( |10|) and the B(rji) are given by (i 
1, 2 , 3, a ik = en - a k ai = b, a 2 = a) 



k 2 - I k 2 - I 

i?(r/ t ) = a 21 ^^+a 12 ^^ 



(176) 



7?i - a r/j - 6 

The wave-functions in parabolic coordinates are given by (£ is the parabolic separation constant) 



n ' c ' pl0 '^'^ ™ w 27rvpr 2 (i + A) 

X^-iC/2p,A/ 2 (-iK 2 )M-iC/2p,A/2(-iCT 2 ) 



(177) 
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(ii) k\ — \ = 0: We have in circular parabolic coordinates 



z(t")=z'' 



V (t")=v" €(*")=e" 

J V V (t) J vat)(f + V 2 ) J T>z(t) 

'Mz'z" 



x exp 



(« 2 + r/ 2 )(e 2 + ? ) 2 ) + i 2 )-^('^i + /C32 ' 



2M\ 



alt} (178) 



M 



exp 



'< ihT 

r,(s")=v" 



V(0)= V ' 



x exp < - 

« Jo 



f ((e 2 + ^ 2 ) + ^(e 2 + r/ 2 )-n 2A;i 



2 1 



1 1 



2M V? 2 ?? 2 

/"OG /■ /'CO 

= / dp z d( d P e-^ h % tX „{z , ,?,T/)% z ,U* , ',? , ,rf') , 

J « IR, 1/ 

with the wave-functions given by (Q is the parabolic separation constant) 



ds 



E 



(179) 
(180) 



x^-iC/2 P ,±fc 1 /2(-^ 2 )M_ iC/2p>±fel/2 (-W) , (181) 

(182) 



3.2.2. We consider the potential (x 3 = z € IR, &i j2 > 0) 

We write down the corresponding path integral formulations 

K^\x",x';T) 
Cartesian Coordinates: 



x(t")=S" t 

J Vx(t) exp J i J y(^ 2 - w 2 (x 2 + x 2 + 4x 2 )) - 

x(t')=x' ^ 



/i 2 /A: 2 



2 1 j.2 1 
K 2 — 4 



2M I x 2 



+ 



(183) 



dt 



(184) 



Parabolic Coordinates (A = 2n ± fci ± fc 2 + 1): 

x exp [A jT [^((e 2 + ^xe 2 + » 2 ) + </> 2 - *w + (e 2 - ^ 2 ) 2 )) 



n 2 



ft l 4 _|_ ft 2 4 



1 



2M{; 2 7] 2 \ cos 2 sin 



(it 



(185) 
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n=0 

tj(s")=»?" €(•")=€" 



7?(0)=r?' 



c(o)=e' 



x exp 



ft 



r" f f + ^ - ^ 6 + ^) + ^ 2 + ??2) - %2 ^ur (f + 1 



2M Vr ?r 



(186) 



Circular Polar Coordinates: 



P (t")=p" 



<P(t")=<P" 



z(t")=z" 



P Vp{t) / Vcj>{t) / Vz{t) 



P {t')=p> 



xexp - 



M 



^( P 2 + P 2 ^ 2 + i 2 - w V + 4z 2 ; 



h 2 f kf j k\ j 1 



2M p 2 \ cos 2 cj) sin 



(it I 
(187) 



Circular Elliptic Coordinates: 



r,(t")=r?" 



*(*")=*" 



/" Z>f(t) /" Pr/(t)d 2 (sinh 2 £ + sin 2 r/) y 2?z(t) 
expj-y^ — d 2 (jsinh 2 £ + sin 2 r7)(£ 2 + 



x 



z(t')=z' 

r) 2 ) — w 2 (cosh 2 £ cos 2 ij + sinh 2 £ sin 2 rj) + 4z 2 



p _ l 



p _ 1 
I , K 2 4 



2Md 2 I cosh 2 £ cos 2 r/ sinh 2 £ sin 2 77 



(188) 



and <3?! a '^ denote the wave-functions (88). K^ V2 \T) is explicitly soluble in cartesian and circular 
polar coordinates by means of the harmonic oscillator in z, and the radial harmonic oscillator 
in x\,X2, respectively (note the similarity to the Holt potential) 

K^\x",x';T) 
Cartesian Coordinates: 



Mlo 



ih sin 2ujT 



cxp 



Mlo 



ih sin 2toT 



(z' 2 + z" 2 ) cos 2loT - 2z'z'' 
2ih 



Mlo 



i?isin loT 



2 2 



II V x i x i ex P 



/ Mojx\x'1 
liftsincjT 



— 6 N ^ $n z ,ni,n 2 { z > ^1) x 2)^n z ,n 1 ,n 2 ( z i x l, x 2) 

n z ,ni,ri2=0 

Circular Polar Coordinates (A = 2n ± &i ± &;2 + 1): 



(189) 
(190) 



ift sin 2loT 
Mlo 
ih sin wT 



cxp 



i/i sin 2a; T 



(z /2 + z" 2 ) cos 2wT - 2z'z h 



]T $i ±fc2 ' ±fcl )(0 , )^ ±fc2 ' ±fcl) (^") exp 



n=0 



— (p + p )cot^T 



\i?isin wT 



(191) 
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£ e -iE N T/hy r 
n z ,n,m=0 



l (^V)*n,,n I m(A0V') , 



(192) 



where the wave-functions and the energy-levels are given by (N is the principal quantum number) 

Cartesian Coordinates: 

2Mu Jj (Muj\ ±k ^ 2 I 



(xi,x 2 ,z) 



l/2±fci 

x a; • exp 



n,;! 



r(ni ± h + 1) 



2ft V 



2Mu 1 
\| V vrft 2 n *n z \ 



exp 



Mtti 



ft 



'2Mw 



ft 



.Btv = M-ZV ± fci ± k 2 + |) , JV = 2(m + n 2 ) + , 
Circular Polar Coordinates (A = 2n ± fci ± &2 + 1) : 



(193) 
(194) 



= *£ ±A!3 ' ±fcl >(0) 

\2Mu fMoj 9 \ x / 2 f 



'2Mw 1 



\ V vrft 2 n *nJ 



exp 



77? ! 



r(m + A + 1) 



exp 



ft 

^2ft~ 



£ 2 



'2Mw 



E N = Tiuj(N ±k 1 ±k 2 + j) , N = 2(n + m) + n z , 



(195) 
(196) 



and the ^n^ 2 '^ ((f)) denote the wave-functions fl88|). 

5.5.5. We consider the three-dimensional potential (ki^ > 0) 

.2 / U 2 l ^.2 _ l 



Vdx] 



a 



^x 2 + y 2 + z 2 2M\ x 



(197) 



This potential can be interpreted as a maximally super-integrable generalization of the Coulomb 
potential. Note furthermore that the pure Coulomb potential has the same separating coordinate 
systems as VJj(x) and therefore the additional centrifugal terms do not spoil this property. The 
pure Coulomb path integral was discussed by many authors, e.g. Castrigiano and Stark [15], 
Chetouani and Hammann |lj| [2^], Duru and Grosche (5f|, Kleinert |3^, [[!]], Inomata [3£], 
Kleinert |$4j], Pak and Sokmen |109| ], Steiner [ 118 1 , and Storchak [f!22[| . We formulate the path 
integrals in which this potential is separable 



K^{x",x';T) 

x(t")=x" 

Vxit) exp < — 
I ft 

x{t')=x' K 

Conical Coordinates: 

r{t")=r" a(t")=a" 



M ■ 2 a h 
x -I > 

o 9 A.f ^ 



2 2 U 2 _ 1 

"i 4 



r 2M frf x 2 



(198) 



/3(i")=/3" 



/" r 2 £>r(t) /" Da(t) y £>/?(i)(/c 2 cn 2 a + k' 2 cn 2 (3) 



r(t')=r' 

x exp 



a(t')=a' 



if 1 M / 

ft if TV 



r 2 + r 2 (A; 2 cn 2 a + k' 2 cn 2 (3)(a 2 + /3 2 ) ) + - 



a 
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k 2 - I k 2 - I 

rc l 4 _|_ ^2 4 



2Ml S n 2 adn 2 /3 cn 2 acn 2 /3 



(199) 



Spherical Coordinates: 



r(t")=r" 

= J r 2 Vr(t) J sm6V9(t) J V<j){t) 

r(t')=r 



6(t")=e" 



4>(t")= 



o(t')- 



x cxp • 



M 



4>(t')=4>' 
r 2 + r 2 9 2 + r 2 sin 2 6<j) 2 ) 



a 
+- 



n 2 



k 2 - 1 t 2 - 1 r 

K l 4 , K 2 4 1 



Parabolic Coordinates: 

r,(t"W €(*")=€" 



expU^" [f ((e 2 + ^ 2 )(e 2 + ^ 2 ) + e 



j?(t')=^' 

x 



r 2Mr 2 sin 2 # I cos 2 . 
4>(t")=<fi" 

I 

4>(t')=<t>> 



+ 



sin 



dt 



^ 2 <P 2 



+- 



2a 



£2 + ^2 2Mi 2 T) 2 \ COS 

Prolate Spheroidal II Coordinates (Ai = 2n ± k± ± + 1): 



- ' ■ sin 2 (6 4 





i/(t»)= 


/ 2? M (t) 


/ 


M (t')=M' 


i/(f)= 



XV(i)(i 3 (sinh 2 + sin^ z/) sin z/ sinh// 



<J>(t>)=<p> 



exp | ^ ^ ~^"^ 2 (( srn h 2 + sin 2 v)(fi 2 + z> 2 ) + sinh 2 n sin 2 z/</> 2 ^ 



'h 2 -l k 2 - I 

K l 4 + K 2 4 



(200) 



(201) 



1 



d(cosh n + cos z^) 2Mii 2 sinh 2 /j, sin 2 ^ y cos 2 4> sin 2 4 



(d 2 sinh // sinh M " sin z/ sin i/") _1/2 £ <^ ±fc2 ' ±fcl) (^')4 ±fc2 ' ±A;i) (<£') 



2vr?i 



-iET/h 



dtj 

(202) 

/•oo 

/ ds" 
Jo 



*»(«")=**" 



v(s")=i>" 



Vix{s) 



Vv(s) exp ■ 



Mo)=/V 



M 



(/i 2 + z> 2 ) + Ed 2 (cosh 1 n - cos z v) 



v{0)=v' 

+ad(cosh n — cos v) I 



2MUinh 2 / u sin 2 f 



(is 



(203) 



We now obtain (Ai = 2n ± &i ± /C2 + 1, A2 = m + Ai + ^, k = a^/—M/2E/h, only the radial 
Green function can be explicitly evaluated) 

4 / dTe iT ^( y3 )(f",x';T) 
n Jo 

Spherical Coordinates: 
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m=0 



~nh*\k E »~ E J* P h 2 p 2 /2M - E 



(204) 
(205) 



Parabolic Coordinates (lo = \J —2E jM\. 

= (^) E*n ±fc2 ' ±fcl) (^')4 ±fca ' ±fc1 '' ' " 



7 n=0 

^ / Mloti't)" \ / Mw^' 
1 V ifi, sin w.s" / 1 \ i??, sin w.< 



sin ws" 
Mto 



00 



n=0 Lni,n2=0 



\ in sin cos" ) V ift sin ojs" / 6X ^ _ 
— ' En — E 



2ias"/h 

2 



(206) 



+ 



1 U 2 p 2 /2M-E 



Jo Jtr 



The discrete state wave-functions in the spherical coordinates 



(207) 



are 



*WM,r) = 4 ±fe2 ' ±fcl) (0)i/ (m + Ai + ^^^^(cos f ) 



x 



(n + Ai + i) 2 
x exp 



m! 
2/! 



a 3 (Z + A 2 + i)r(/ + 2A 2 + 1) J \a(l + A 2 + ±) 



1/2 



2r \ A2 



a(Z + A 2 + 



IT H 



(2A 2 ) 



2r 



-E^ = — 



Ma 2 



,a{l + \ 2 + \)J ' 
-5— - N = l + m + 2n±k 1 ±k 2 + l ■ 

h A N 2 

The continuous wave-functions are 

*WM,r) = (m + \i + \) T{m + ^^ PxX( cos0 ) 

T(i + A 2 -i/ap) / 7T \ . . 

X v^rr(2A, + l) eXP (^) M ^ ( - 2lPr) - 

The discrete state wave-functions in parabolic coordinates are 



(208) 
(209) 



(210) 



2ni!n 2 ! 



*w«>M = ^ fci) ^[^- r(n 1 + A 1 ;iir ( ; 2 + A 1 + i) 

V (aN) 2 J eXP \ 2aN 

Ma 2 

En = --0-^7, N = m + n 2 + 2n±k 1 ±k 2 + l 
The continuous wave-functions are 



1/2 



(212) 



,m, m, ^ r^ + ^ i/a + c r ^ + ^ l/a-C ,„ 

B * (WMJ n W 2vr v ^^r(l + Ai)r(l + A 1 ) 

xM_ i(1/a+0/2piAl/2 (-ip^ 2 )M_ i(1/a _ c)/2piAl/2 (-ip?7 2 ) . (213) 
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For the expansion into the wave-functions one uses ( |135|) and performs the variable substitution 
(p£,Prj) >— > [(V° + O/^P) (l/o — C)/2p], a = h 2 /Ma. ( is the parabolic separation constant. 



We consider the potential (£^2,3 > 0) 



Va(x) 



kix 



2M \y 2 v / x 2 + y 2 y 



k 2 - I k 2 - I ' 
+ ^2 4 _j_ ^3 4 



2" 



(214) 



Note that for /cf - \ = we obtain a trivial extension of the a = case of the two-dimensional 
maximally super-integrable potential V^{x). We formulate the path integral in spherical, circular 
polar, circular parabolic and circular elliptic II coordinates, respectively 

K^\x",x';T) 



x(t")=x" 

Vx{t) exp < — 

x(t')=x' ^ 

Spherical Coordinates: 



k\x 



M ^2 ft 

Y X ~ 2M Wz 2 + V 2 



2 _ 1 1.2 _ 1 



(it } (215) 



r(t")=r" 8(t")=e" 4>(t")=<t>" 

„2 



r^Dr(t) 

r(t')=r' 0(f) 



sin 6V9(t) 



0(t')=0' 



xexp|- 
ft 2 



y(r 2 +r 2 2 + r 2 sin 2 6»0 2 ) 



1 /* H~~ 



2Mr 2 lsin 2 0V 4sin 2 ; 
Circular Polar Coordinates: 

P {t")= P " <P(t")=<P" z(t")=z" 



+ 



k\ k\ \ 1 



4 cos 2 ^ 



T + 



^1 I _ 1 

cos 2 9 4 



eft 



Vz(t) 



p(t')= P ' 



xexp - 



<t>(t')=4>' 



M ,-2 



z(t')=z> 



212 . 2) 



2M 



1 / 1,2 _ h.2 _ 1 p 1 p _ 1 



P V 4 COS 2 



+ 



4 sin 2 I 



1 



+ 



1.2 1 

rc 3 4 



Circular Elliptic II Coordinates: 



£(*")= 



r?(t")=r?" 



z(t")=«" 



/" X>f(t) / P? ? (t)ci 2 (sinh 2 ^ + sin 2 ? ? ) / 2?z(t) 



€(f)=C' 
x exp 



•7(f) 






' M / 




y' 



z(*')=Z' 

sinh 2 £ 2 + sin 2 77) (£ 2 + ?) 2 ) + i 2 ) + 



n 



ci(cosh£ + cos 77) 



/V 2 



1 



2 cosh £ cos jy + l 2 

1 t — 7 2 



1 



2M I d 2 sinh 2 £ sin 2 77 V 1 cosh £ + cos 77 
Circular Parabolic Coordinates: 

z(i")=z 



+ 



A- 2 - i' 

ft 3 4 



(it 



»7(t"W 

/ 

7j(i')=r/' 



e(*")=e 

| Pr?(t) | pe(t)(e 2 + r? 2 ) | 2?if(i) 

e(f)=f 









{{ 


r 





z(t')=«' 
2\/^2 , -2\ . -2 



(216) 



(217) 



(218) 
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Ti 1 



+ 



We obtain for the path integral solutions (A± 
2m + Ai ± k 3 + 1) 



+ 



dt 



(219) 



n + (A+ + A_ + l)/2, A 2 
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Spherical Coordinates: 

= -(rV , sm^sinO -1/2 E^ A+ ' A - ) f-)*l A+,A - ) 



M 
x -ir^; exp 

ihT y 



n=0 

[M y 2 +r» 2 ) 



2hT 



m=0 



ihT 



00 fOO 



(220) 
(221) 



n,m=0 

Circular Polar Coordinates 



z' z" exp 



,2 1,2 ,2 „2 
— {z +z +p +p 



7 ±fc 3 

Mp'p" 
ihT 



V ihT 



n=0 

00 />oo /-00 

n=0 io Jo 



(222) 
(223) 



Circular Parabolic Coordinates: 



M n-r, 

— — V z'z" exp 
ihT 



^L(z' 2 + z" 2 ) 
2hT y ' 



Mz'z'' 
ihT 



x / d( 



- ^ |r(^ + |)| 2 |r(^ - |)| V/ Q p 
p 4vrVe , ewr 2 (i + A_)r 2 (i + A+) 

t«2 



-ihp 2 T/2M 



xA^_ iC/2pIA+/2 (-^ , ")M i^/2plA+/2 (ipe , )M iC/2piA _ /2 (-W z )M_ iC/2piA _ /2 (W z ) (224) 



R JO 



dp, / dc/ dpe- iB — r /^; jiPiC (^,e / ,7 7 / )^ >P ,c(^,e",V / ) 



(225) 



(C is the parabolic separation constant) where the wave-functions and the energy-spectra are 
given by 



Spherical Coordinates: 



Ep= 2M P ' 

Circular Polar Coordinates: 



(2rsin^)- 1 /2^ +) A_)/|j $ (A 1> ± fc3 ) Wv ^j A2(pr) 



E 



(226) 
(227) 

(228) 
(229) 



Circular Parabolic Coordinates: 



^Pz,P,d z >^v) = VpTz J±k 3 (p z z) 



TV 1+A + 1 'C \-p r 1+A - _ Kv/ 2 * 

1 V 2 ~r 2p' L V 2 2rJ e 



2p/ 



27Tv^r(i + A_)r(i + A + ) 

x^-iC/2p,A + /2(-ipC 2 )M_ if/2piA _ /2 (-ip? ? 2 ) 



(230) 
(231) 



In the path integral evaluation subsequently [55| the path integral solution for the Poschl- Teller 
potential and the radial harmonic oscillator have been used. The spectral expansion in the 
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z- and p-dependent kernels is done by means of (analytical continuation [113] of) the Weber 
formula 



re iar I v (—\ar)I v (—ibr)dr 



_j_(a 2 +b 2 )/4ai j ( a b 



2a 



\2ai 



(232) 



and we obtain the path integral identity (we use the functional weight formulation) 

r(t")=r" 



r 2 dt 



rrr, M 

yr'r"— — exp 



iM , , 2 



i //2n 

r + r 



r(t')=r' 



(Mr'r"\ , . 

V^V 7 / pdpe- ihpT/2M J x (pr')J x (pr") . (234) 
Jo 



Applying (|234| ) gives in an obvious way the spectral expansion of the z- and p-dependent kernel. 
For the extraction of the continuous states in parabolic one applies ( p.35| ) by introducing a 
Coulomb coupling a and sets a = oo, i.e. infinite Bohr radius, in the in the final formula?. 

In circular elliptic II coordinates the z-path integration is separated in the same way as for 
the circular parabolic case, and for the remaining path integrations we obtain the identity 



K^(x",x';T) 



M 



VTz 7 ' 



exp 



iM <z> 2 + z" 2 ) 



2hT 



l ±k 3 



Mz'z'' 



ihT 



SR. 



dE 
2vrft ( 



t{s")=Z" 



v(s")=v" 



xexp - 



AET/h 



M 



ds h 



Vt](s) 



e(o)=e 



J7(0)=7?' 



y (£ 2 + V 2 ) + Ed 2 (cosh 2 £ - cos 2 rj) 

ft 2 f {k 2 - I) + kfcosT] (k\ - \) + A; 2 cosh £ 



2M 



sin 2 r\ 



+ 



sinh 2 £ 



ds\ . (235) 



For the spherical coordinates the expansion into radial wave- functions, and for the circular 
parabolic coordinates the expansion of the ^-dependent kernel are performed by means of the 
Hille-Hardy formula fl93|). 



3.2.5. We consider the potential (kip > 0) 



V 5 (x) 



h 2 



k\x 



+ 



k 2 - i' 

"-2 4 



2M \y 2 v / x 2 + y 2 y 



k*z . 



(236) 



Note that for k% = we obtain a trivial extension of the a = two-dimensional maximally 
super-integrable potential Vs(x). We obtain the path integral formulations 

K<y*\a",x';T) 

Circular Polar Coordinates: 



p{t")= P " 



4>{t")=4>" 



z(t")=z" 



P Vp(t) 



V<j>(t) 



Vz(t) 



p(t')= P > 



xexp - 



<t>{t>)=<t>> 



z(t')=z' 



— (p + p (j) +Z ) 



kj,Z 



' k 2 - k 2 - I k 2 + k 2 - 1 



2Mp 2 



4 cos 2 | 



+ 



4 sin 2 | 



dtj 
(237) 
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Circular Elliptic II Coordinates: 



/ 



I7(t")=»j" 

/" 2?^(t) J Vri(t)d 2 {smh 2 ^ + sin 2 rj) / X>*(t) 
»?(t')=J7' 



z(t")=z h 



z(t')=z> 

xexp I ^y y((sinh 2 ^ 2 + sin 2 ?7)(^ 2 + ?) 2 ) + i 2 ) - k 3 z 



l 



2 cosh £ cos 77 + I 



2M d 2 sinh 2 £ sin 2 77 ^ 1 cosh £ + cos 77 
Circular Parabolic Coordinates: 

z(t")=z r 



+ fc 2 — I 



(it 



(238) 



/ 

ri(t')=V' 



xexp - 



«(*")=? 

/ 

?(i')=e 

M 



| 2?r?(i) | P£(t)(£ 2 + r7 2 ) | 2>*(t) 



*(#)=«' 



^((e 2 +r7 2 )(£ 2 +r7 2 ) + i 2 )-- 77 (^ 



Ti z 



v 2 -e 



+ 



ft 2 4 



2M\ i e 2 77 2 (£ 2 + r7 2 ) £V 



(239) 



Parabolic Coordinates (A 2 = /c 2 + A; 2 ) 
r?(t")=rj" £(*")=£" 



V(t')=r,' 
X 



0(i")=0" 

| 2?r?(i) I P£(t)(£ 2 + r7 2 )£r7 | 2ty(i) 

0(i')=</»' 



exp jijf y((e 2 +r7 2 )(£ 2 +r7 2 )+e 



2 212 

77 



f*l+*2 4 _ fc 2_ I] _ ^3(^2 ,^2) 



2Mf V ^ sin 2 1 4 

(rtW)- 1/a f> + A + i) ^ + 2A + 1 ) p- + A Jcos0O p- + A n( 



(240) 



COS ( 



n=0 



»7(«"W 



IR 



2tt}V 



€(•")= 



e -iET/n / ds n 



V (0)= V > 



Vn(s) J V£(s 

«(o)=e 



xexp - 



n 2 ( \ 2 -k 2 -\ | A 2 -fc 2 -f 



2M 



ri- 



ds 



(241) 



The path integral in parabolic coordinates cannot be done due to the £ 4 , 77 4 dependence. Actually 
V5 represents a Stark effect. In the cases of the circular polar and circular parabolic coordinates 
the corresponding Green function can be evaluated by means of the path integral solution for 
the linear potential for the z-path integration, and the radial harmonic oscillator in the 
(r, (p), and in the £- and r7-path integration [A 2 ^ = k 2 ± k 2 , Ai = n + (A + + A_ + l)/2] 

K^\x",x'-T) 

Circular Polar Coordinates: 



( m y/ 2 

{^mr) exp 



i IM , 2 k 3 T „ s k 2 T^ 



2AM 
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M 
2\KT 



n=0 



^ $ (A +1 A_)/^_U(A + ,A_) 



exp 



" (/ 2 + / 2 ) 



2\%T 



hi 



Mp'p" 
ihT 



00 poo f 

n=0 Jo Jm 



(242) 
(243) 



Circular Parabolic Coordinates: 



( m y/ 2 
l^rj exp 



ft 2T U 



24M 



x / dC 



- d P \rC-^ + |)l 2 |r(^ - 1)| 2 ^ /a \ _ mv , TI2M 



In Jo P 4vrV^'Wr 2 (l + A_)r 2 (l + A+) 
xM_ K/2pM/2 (-i P e 2 )M K/2p!X+/ ^^ (244) 

/■ /'(DC- r 

= d( dp dE z e-^ T /^ PtC (z\eW)^E z , P x^'\ey) , (245) 

(C is the parabolic separation constant) where the wave-functions and the energy-spectra are 
given by 



Circular Polar Coordinates 

, 2M y/ 3 



h 2 Vk 



Ai 



h 2 



Ee v = E z H p z , 

,p 2M 

Circular Parabolic Coordinates: 



2M X 1 / 3 



h 2 Vk 



Ai 



E z \ [2Mk\ 1 / 3 



E 7 \ /2Mk\ 1 ^ 



2-V2$(A + ,A_)(|j v ^ (7Ai(pp) )(246) 

(247) 



A; 



r /i+A+ . iC frv 1+A - — Kv 7r / 2tt P 

1 V 2 ' 2p/ x V 2 2r>>^ 



2p> 



Ee d = E z -\ p 

z ' p 2M 



27rv^?^(l + A_)^(l + A^ 

h 2 



-M_ K/2pM/2 (-ipe)M K/2ptX _ /2 {-W) , (248) 



(249) 



C is the parabolic separation constant. 

The expansion into the wave-function in the z coordinates is done according to ( ]3"q] , Ai(z) 
are Airy functions) 



z(t")=z" 

j T>z(t) exp 

z(t')=z' 

( m y/2 

= \2^T) 6XP 



m 



z 2 — kz\dt 



h J t > V 2 



n \2T 

2M \ 2 / 3 



k 2 T 3 \ 

2 v " ' " ' 2AM ) 

. E\/2Mk\ 1 / 3 , 
11 ' Ai 



k)\ h 2 J 



E\ /2M£;\ 1/3 



kj 



(250) 
(251) 



The expansion into the radial wave-functions for the circular polar coordinates is performed by 
the Weber formula (234). 



3.3. Three-Dimensional Minimally Super-Integrable Smorodinsky-Winternitz Po- 
tentials. 



We discuss in this subsection the eight three-dimensional minimally super-integrable potentials. 
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The are characterized by having four functionally independent integrals of motion. In the sequel 
x denotes a three-dimensional coordinate: x = (x,y,z) = (xi,X2,xs). 

In table 4 we list the three-dimensional minimally super-integrable Smorodinsky-Winternitz 
potentials together with the separating coordinate systems (where the italized coordinates sys- 
tems were not mentioned in |35|]. The cases where an explicit path integration is possible are 
underlined . 

Table 4: The three-dimensional minimally super-integrable potentials 



Potential V(x, y, z) 


Coordinate System 






Spherical 
Conical 


rr M 2, 2 2n & f k l ~ \ 


+ 2^ 


+ F(z) 


Cartesian 

Circular Polar 
Circular Elliptic 


V 3 = A *^(4x* + y 2 ) + f k \hF(z) 
2 2M y A 


Cartesian 
Circular Parabolic 


Vi- a n2 ( h * x 

yjx 2 + y 2 2M\y 2 ^J x 2 + y 


k 2 — 

2 y2 




Circular Polar 

Circular Elliptic II 
Circular Parabolic 


_ j. M 2 2 2 2^ K 2 f k 3-J F{y/ X )\ 


Spherical 

Circular Polar 
Oblate Spheroidal 
Prolate Spheroidal 


M 2/2 2 , 2n ^ F(y/x) 

* = 2^ X +y + ^ + 2Mx 2 ly 2 


Circular Polar 
Parabolic 


v - a i n * ( k2lZ 1 f( 

7 " r 2M(x 2 + y 2 ) \ r \ 




Spherical 
Parabolic 

Prolate Spheroidal II 


a p2 f 4> 4> s 
% = -- + y-(ftco B - + fc B m- 




Mutually Circular 
Parabolic 



3.3.1. We consider the potential (r = &i,2,3 > 0) 

f, 2 fh 2 - i h 2 — - h 2 - I\ 
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We formulate the path integral in the separating coordinate systems and obtain 
K^\x",x';T) 



x(t")=x" 

T>x[t) exp < - 

x{t')=x' ^ 

Spherical Coordinates: 



M J2 . . 
— x — Fir) , 
2 v ' 2M\ x 2 



21 k 2 1 u2 1 



n / ft l 4 i ^2 4 i 3 4 



+ 



+ 



(it ^ (253) 



r(t")=r 

/ 

r(i')=r' 



6(t")-- 



<P(t")=<P" 



J r 2 Vr(t) J sm8V9(t) J V<f>(t) 



xexp - 



6(t')=8' 



M 



— ir 2 + r 2 6 2 +r 2 sin 2 OS 2 ) 



F(r) 



l /k( 



2Mr 2 \ sin 2 6 V cos 



2 I 1.2 _ l 

_4 , ^2 4 

, T" . o , 



1\ A; 



sin 



+ 



2 1 

3 ~ 4 



1 



cos 2 9 4 



(254) 



Conical Coordinates: 



r(t")=r" 

= J r 2 Vr(t) J Va(t) J V(3{t)(k 2 cn 2 a + k' 2 cn 2 (3) 

r(t')=r' 



a(t")=a" 

I 

a(t')=a' 



0{t")=(3" 



— [r 2 + r 2 (fc 2 cn 2 a + k' 2 cn 2 [5){a 2 + >r ) ) - /-•'( , ) 



k 2 - I 

K 2 4 



2Mr 2 ysn 2 adn 2 /3 cn 2 acn 2 /3 dn 2 a;sn 2 /3 



^3 4 



(it 



(255) 



For this potential the angular path integrations can be easily done and one is left with a radial 
path integral pOfl, i.e. (Ai = 2n ± k\ ± &2 + 1, A2 = 2m + Ai ± ^3 + 1) 



K^\x\x'-T) 

00 00 
= (r'V 2 sin 0' sin #"r 1/2 E ^ k ^\<p")^ k ^ ±k ^ ((/>') £ ^'^(fl")^ 1 '^ 3 ^') 

r(t")=r" 



n=0 



m=0 



X 

r(t')=r" 



Vr(t) exp 



i /•*" /M. 2 n 2 A 2 -f 
n A' \ 2 r U 2M r 2 



(256) 



(r'V 2 sin & sin 0~ 1/2 E ^ ±fe2 ' ±fel )(^")^ ±fc2 ' ±fcl) (</>') E ^'^(O^ 1 '^ 3 ^') 

n=0 m=0 

r(t")=r" 



ji\ 2 [r 2 ]£>r(t) exp 



r(t')=r" 



00 „ 

Y, y^A I .e- i ^/^ Lin)m (^,0V , )^A r ,n im (^,0 ,, ,r'') . 



(257) 



(258) 



m,n=0 



Here use has been made of the functional weight formulation in path integrals according to 
H The wave-functions have the form 



(r 2 sin 0)-V 2 $(±fc2,±fci) ^^(Ai,±fc 3 ) ^ 



(259) 



where $A r ( r ) are the radial wave- functions. Of special interest are the cases F(r) = mw 2 r 2 /2 
and F(r) = —ot/r, c.f. 3.2.1 and 3.2.3, respectively. 
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3.3.2. The next three potentials represent a simple z-dependence extension of the two-dimensional 
potentials V1-V3. We consider the potential (k± 2 > 0) 



m n 2 ( k 2 - I h 2 - i' 



+ 



+ F(z) . 



(260) 



For the case of F(z) = 2Mu?z 2 there is an additional coordinate system, the parabolic coor- 
dinates, which separate V2, c.f. the three-dimensional maximally super- integrable potential V2. 
Here now we have 

K^(x",x';T) 
Cartesian Coordinates: 



S(t")=S" 



x{t')=x> K 

Circular Polar Coordinates: 



M C-* 2 2/2,2 

— [x -u [x +y 



n 2 2 k 2 — - 

' v * 4 - F(z) 



E 



2M^[ 



dt 



(261) 



z(t")=z' 



p(t")=p" 



Vz{t) / pPp(i) 



z(t')=z' 



p(t>)=p> 



xexp - 



/" -2 , 2i2 2 2\ 
— (p +p - W p ) - 



?j 2 / fc 2 — ^ fc| — j 1 



2Mp 2 \ cos 2 <fi sin 



^ (262) 



Circular Elliptic Coordinates: 

z (t")=z" £(*")=£" 



y D«(t) J V£(t) J Pr/(t)d 2 (sinh 2 ^ + sin 2 r/) 

exp < ^ jf — (^d 2 (sinh 2 £ + sin 2 r?)(£ 2 + ?) 2 ) + i 2 — u; 2 d 2 (cosh 2 £ cos 2 77 + sinh 2 £ sin 



p _ 1 L.2 _ 1 

K l 4 _|_ ft 2 4 



2Md 2 I cosh 2 £ cos 2 r\ sinh 2 £ sin 2 77 



r/j 
(263) 



Except for the z-path integration, we can perform in the case of cartesian and circular polar 
polar the path integration explicitly yielding 



K<y»\£",x';T) 
Cartesian Coordinates: 



Muj 



ih sin ojT 



3 2 



n v x i x i ex p 



i=i 



Mu (x? + xf)cotuT 



z(t")=z" 



Vz(t) exp 



z(t')=z' 



2ih 



hi 



( Mux\x" 
Viftsin uiT 



00 „ 



(264) 



(265) 



ni,7i2=0 

Circular Polar Coordinates (Ai = 2n ± &i ± &2 + 1)): 



i?isincjT 

xK F (z",z';T) 



g $ (±fc 2 ,±fc 1 )(0')$(±fc2,±fcl) ( ^/ ) exp 



n=0 



M(J /2 „2x , „ 

— (p +p )cot^T 



/ Mwp'p" 
Kl \ihsmujT 



(266) 



46 



oo - 
n,m=0 

and in the sequel we denote the remaining z-path integral by 

K F (z",z>;T) = J dE Xz e- iE ^ n n z (z')^ Xz (z") . 
The wave-functions in (xi,xz) and (p, (j>) are the same as in 3.1.1. 
5.5.5. We consider the potential (k 2 > 0) 

M 



(267) 



(268) 



+ F(z) . 



(269) 



Due to the similarity of the corresponding two-dimensional case in 3.1.2 we just give the path 
integral formulations which have the form 

K<y»\x",£';T) 
Cartesian Coordinates: 



x(t")=S" 



J Vx(t)ewU£ y(f 2 -^ 2 (4x 2 + y 2 ))-F(z) 



^ 2 ^2 — I 



2M y 2 

T ' h * n ,iwy)* n ,iW>i/) I dE ^ e- iEx * T/n n z (z')*x z (z") 



x(t')=x' 
oo 

n,l=0 

Circular Parabolic Coordinates: 



dt 



(270) 
(271) 



v(t")=v" 

I 

r,(t')= V ' 



x exp 



M 



k 



e + v 2 )(e + v 2 ) - ^ae - v 2 r - ev) )-^(e- v 2 ) - n 



2 fc 2 I 

2Mr) 2 i 2 



xK F (z",z';T) , 

with the same Feynman kernel and wave-functions in (x,y) as in 3.1.2. 
3.3.4- We consider the potential (k\^ > 0) 

h 2 



(272) 



V 4 (x) 



a 



+ 



\/x 2 + y 2 2M \ y 2 \Jx 2 + y 



k\x 



+ 



^2 4 

2 

y 



(273) 



We formulate the path integral in circular polar and circular parabolic coordinates, respectively 
K^I(i",x';T) 



x(t")=x" 

- 1 

x(t')=x' 

Circular Polar Coordinates: 



a 



h 2 



k\x 



M % 2 

T X + ^Jx 2 + y 2 2M\y 2 ^x 2 +y 2 ' y 2 



f ^M^ - F(z) 



(274) 



P (t")=p" 

I 

P (t>)=p> 



J pVpit) J V4>{t) 

4>{t>) 
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x cxp 



h Ji 

Jl J 



M.f-2 2 A?\ a V [ k\-k 2 -\ k\ + k\-\ 
2 [p+P ^ )+ p 2MA 4cos 2 ^ + 4sin 2 4 4 



1 



dt 



xK F {z",z;T) 
Circular Elliptic II Coordinates: 



(275) 



an-- 

I 



V (t")= v " 

J V£(t) J Vrj(t)d 2 (sinh 2 £ + sin 2 rj) J Vz(t) 
ri(t')=n' 



z(t")=z" 



x exp 



z(t')=z' 

^(( s inh 2 ^ 2 + sin 2 r ? )(C 2 + ? ) 2 )+i 2 ) + 



a 



d(cosh£ + cost/) 



n 2 



i 



2 cosh£cos?7 + l 2 l 

rCi " " — (~~ Kn 



2M d 2 sinh 2 f sin 2 r, \ 1 cosh £ + cost] 2 2 
Circular Parabolic Coordinates: 

z(t")=z" 



F(z) 



dt } (276) 



v(t")=v" 

I 

T](t')=rr 



£(*")= 



| 2?7/(t) | P^X^ + r/ 2 ) | X>*(t) 



» p I s jf If (« 2 + " 2 »« 2 + " 2 » + - i2 ) + wh - F( --> 



rf-e 



+ 



k 2 - I' 

K 2 4 



2Ml ^V(e 2 + r/ 2 ) £V 



(277) 



The circular polar an circular parabolic path integral formulations can be explicitly evaluated 
in terms of the corresponding Green function (A± = k 2 ± A; 2 , Ai = n + (A+ + A_ + l)/2, 
k = a v / -M/2E/h) 

KW(Z",x';T) 

Circular Polar Coordinates: 



k f (z", A T){ P ' P "T 1/2 \ E ^ A+ ' A - } (f ) 4^ ( y 



1 /• d£ 



-iET/h 



n=Q 

M r(~ + Ai — k) 



ft 7ia 27r?i 



V 2£ T(l + 2Ai) 
*A 1 ,m(P / ^ , )^A 1> m(/,0") 



W Ml ( ^8ME ^ ) M K)Al ( v/^SM^ ^ ) (278) 



P< 



oo oo 
X E E 



E n — E 



m=—0 I n=0 

Circular Parabolic Coordinates: 
= JdE Xz e- iE ^ h n z (z')^ Xz (z") 



+ I dp 



h 2 p 2 /2M - E 



(279) 



E 

ni,ri2=0 



K — E 



roc r 1 

+ / dp d(- 

Jo Jn 



h 2 p 2 /2M -E 



(280) 



with the same wave-functions in ((f), p) and (£, rj) as in 3.2.3. 



3.3.5. The next three potential involve an arbitrary ^-dependence where F(y/x) = 7 2 > 
describes the well-known double ring-shaped oscillator (c.f. e.g. Carpio-Bernido et al. [O, O], 



Kibler et al. |76|, |77|, [81]], Lutsenko et al. |9T|], and Quesne |116| ). We consider the potential 
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(fc 3 > 0) 



M 



V b (x) = ^ 2 {x 2 + y 2 + z 2 ) + ^- 



h 2 k 



.2 1 

'3 4 



+ 



x 2 + y 2 



(281) 



2 v 7 2M 

We obtain the following path integral formulations in the corresponding separating coordinate 
systems 

K<™(x",£';T) 
Spherical Coordinates: 

4>{t")=4>" 



r(t")=r' 



r(t')=r' 

x exp 



= J r 2 Vr{t) J sm9V9{t) J V(f>(t) 
e(t')=e' 
i /■*" \M / 
hit TV 



r 2 + r 2 # 2 + r 2 sin 2 6>0 2 - w 2 r 2 



h 2 (F(t & nct>)-\ k 2 -\ 1 



2Mr 2 

Circular Polar Coordinates: 

z{t")=z" p(t")=p" <t>{t")=<t> 



sin 



cos 2 6* 4 



(it 
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= / Z>*(t) / PVPQ) J V<f>(t) 
P (t')=p> <t>{t')=<t>' 

M ( -2 i 212 i -2 2/2 , 2\\ ^~ 

Oblate Spheroidal Coordinates: 



z(t')=z' 



xexp - 



fi 2 / F(tan 0) - 



+ 



t.2 _ i ' 

% 4 
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A(t")=A 



p(t")=p" 



A(t')=A' 
X 



P(t')=i ? ' 



o(/") O" 

J T^U{t) j Vv{t)dP (cosh 2 p, — sin 2 v) cosh jl sin P / T^<P(t) 
; • ' o(f') o' 



exp j ^ ^ ^-(?^(cosh 2 |U — sin 2 f)(/i 2 + v 2 ) + cosh 2 /i sin 2 i/0 2 

—to 2 (cos 2 /Z sin 2 P + sinh 2 if cos 2 u)^ 



h 2 /F(tan^)-± fc 2 -i 



2Md 2 \y cosh 2 ii sin 2 f sinh 2 ju cos 



Prolate Spheroidal Coordinates: 



dt 
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M*")=m" 

/ 

p(f)=ji' 



i/(t")=i/" 

J D n(t) j XV(t)d 3 (sinhV + sin 2 v) sinh // sin i/ ^ 

(/(//) i/' ' ' 



xexp - 



^d 2 ((sinh 2 |U + sin 2 ^)(/i 2 + v 2 ) + sinh 2 /i sin 2 
—lo 2 (sinh 2 /i sin 2 v + cosh 2 /x cos 2 i/)^ 



/i 2 / F(tan , 



2Md 2 ^ sinh 2 ll sin 2 v cosh 2 /x cos 2 v 



h 2 - I 

ft 3 4 



(285) 



We consider the spherical and circular polar coordinate cases, the others can be treated in a 
similar way as far as the separation of the (/>-path integration is concerned. We denote the 
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energy-spectrum of the (^-dependent problem by = % 2 \^ i /2M . The most important case 
occurs when F (taxi <fi) = constant. Then the quantum motion in <f> are just circular waves and 
we have in this case ^a,^) = & u (4>) = e 1 ^ j\f 7 2M (v £ 2). This is also the case for the potentials 
Vq and V 7 . We have (A 2 = 2n + ± k 3 + 1) 
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Spherical Coordinates [10]: 

= (r'r" sinO' sme")- 1/2 J dE^tf)*^) 



Mlo 



J2 $ ( n*' ±k3 \e')$l X *' ±k3 \9") - — - exp 



n=0 



iftsin loT 



^(r' 2 +r" 2 ) cot coT 



h 2 T— T 



Mur'r' 1 



ih sinuiT 



r oo 



n,l=0 

Circular Polar Coordinates: 



Mlo 



ih sin loT 



exp 



— — (z +z + p + p ) cot coT 



l ±k 3 



Mujz'z" 
ihsinioT 



- oo 



m,n z =0 



The wave-functions and energy-levels are given by 



Spherical Coordinates: 



hMuj (Mlo^I 2 



h V n J Vr(/ + A 2 + i) 



r A 2 +l/2 



exp 



Mu} . r A L w( Mu} ^ 



2h 



E N = huj(N ± fc 3 + 2) , /V = 2/ + 2n + A^ , 

Circular Polar Coordinates: 

., , _ . ^2M^/Mc^\ ±fc 3/ 2 



h \ h 



■ 



T(n 2 ± fc 3 + 1 
Mw 9 \V 2 



. z l/2±fc 3 exp 



E N = huj(N ±k 3 + 2 



r(m + A + r 



Mlo 

~ ~W 

exp 



L(±*s) 



Mo; 2 N r (A )/Mw 2 



2ft 



TV = 2m + 2n z + A^ . 
For F(tan0) = constant, \ ( j ) = \v\, c.f. also Vq and V7. 
3.3.6. We consider the potential 
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(289) 



(290) 
(291) 



(292) 
(293) 



(294) 



We obtain the following path integral formulations in the corresponding separable coordinate 
systems We write down the corresponding path integral formulations 

K<y*Xx",$';T) 

Circular Polar Coordinates: 



z(t")=z" P (t")=P" w ,- v 

J Vz{t) J P Vp{t) J V<f>(t) 

z(t')=z> p(t')=p> 



4>(t")=4>" 

I 

$(?)=</>> 
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x exp 



v 



M/. 2 2 12 , -2 2/ 2 . A 2\\ % 2 F '( tan <t>) ~ I 

— ( p +p<?> +2; - (p +Az)\-n 



2Mp 2 



dt 
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z(t")=z" p(t")=p" 

x | 2?s(t) | ^[p 2 ]Pp(t)exp 

z(t')=*' P (t')=p> 



JdE x , g e" ii? -™ T /^; 2iA0im (0',p',z')^,A„ m (0 /, ,p",^) 



(296) 
(297) 



Parabolic Coordinates: 



x , e - iET / a / ds" 

/K 27r/i Jo 



ri(s")=v" 



Vr](s) 



x exp 



f ((f + . 2 ) - - 2 (c 6 + . 6 )) + E { e + , 2 ) - (i + 1 



ds 



(298) 



Here we have used the same notation as in the previous case. The wave-functions and the 
energy-spectrum in circular polar coordinates are given by 



! 2Mu 1 



\ V ^ 2 



exp 



h \ ft 

Muj 



-P 



ml 



r(m + A + l) 



'2Mw 



^n z ,m = ftuj(2m + \<f ) + n z + 



3.3.7. We consider (fa > 0) 
F 7 (x) = 



2/ • 



n 



+ 



/r 



, = + F ( y - 

yjx 2 + y 2 + z 2 2M(x 2 + y 2 ) \ ^ x 2 + y 2 + y 2 \x 



, M " 2 \ T {\^)(Mio 2 



(299) 
(300) 

(301) 



For F(y/x) = j 2 > this potential is known as the ring-shaped Hartmann potential (c.f. 



Calogero §§, Carpio-Bernido et al. 0, 0, 0, H, Chetouani at al. Q, Gerry g§, Granovsky 



et al. @], Grosche @, Guha and Mukherjee f^Hartmann @, Kibler et al. @-|7g], [||, 
Lutsenko et al. [90], Vaidya and Boschi-Filho [125], and Zhedanov [131]; compare also the 



connection with a Coulomb plus Aharonov-Bohm solenoid, e.g. Chetouani et al. [18|, Kibler and 



Negadi | |80| |, and Sokmen [117]). We write down the corresponding path integral formulations 
where only the corresponding Green function can be evaluated [A± = A? ± k\ , X± = n + (A+ + 
A_ + l)/2, k = a v / -M/2E/ft] 

| [°° dT^ ET ' n K^\x\x'-T) 
ft Jo 



ft Jo 

Spherical Coordinates: 

r{t")-- 



6(t"> 



W")=4>" 



^J°°dTe iET/h J r 2 Vr{t) J sm9V6(t) J V<j>(t) 



r(t')=r' 



e(t')=8' 
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x exp 



M 



(r 2 + r 2 6 2 + r 2 sin 2 6<j) 2 ) 

a h 2 ( k\ cosO F(tanc/>) - \ 
+7 ~ 2Mr 2 I sin 2 # + sin 2 9 
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(it 
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/-i oo / ft'\ /a" 
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x i / m r(| + Ai-«) 



fry 2£ r(2A 2 + l) 



W KM (V^mE^M KM (V^mE^j (303) 
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rJ^ol^ £]V-£ ilR P 



h 2 p 2 /2M -E 
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Parabolic Coordinates: 



Wo 



dTe iET/k 



-2 , „2 



2?7/(t) 



x exp 
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((e 2 + r/ 2 )(e 2 + ^) + eV^ 2 ) 
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fi2 / 2 _ t2 
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E 
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E 

ni, 7i2=0 
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4 



dt 
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R — R 

Jj ni,n,2 J -' 



h 2 p 2 /2M -E 



+ r dp i d( 

Jo Jwi 

The bound state wave-functions in the polar coordinates are 



n* 1 c(^^V)*n lPl c(^^£ , W , ) 



(306) 



(2sin0)- 1 / 2 ^(^ ) c I >( n A + ,A-)/_ 



Z! 



(/ + A 1 + i) 2 
2r 

a(/ + Ai + \) 



a 3 (/ + Ai + i)r(/ + 2A! + l) 

Ai 



exp - 



1/2 



(2Ai) 



a(/ + Ai + i)/ \a(/ + Ai + i) 



2r 



Rat 



Ma 2 
FiV 2 



iV = / + n + i(l + A_ + A+) 



(307) 
(308) 



The continuous wave-functions are 



r(i + A 2 -i/op) / 7T \ , . 



(309) 



The wave-functions in parabolic coordinates are the same as in 3.2.3, except for the modification 
that A 2 = A 2 . ± k\ instead of \ 2 ± = k 2 ±k 2 . 
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The potential Vj is also separable in prolate spheroidal II coordinates and we have the 
identity 

K^\x",x';T) 

/u(t")=M" v(t")=v" <t>(t")=<j>" 

= J Vfi(t) J P^(t)(i 3 (sinh 2 /i + sin 2 ^)sin^sinh / u J V(j>{t) 

H(t')=ii' v{t')=v> <f>{t')=(t>' 



xexp W 



— d 2 ^(sinh 2 fx + sin 2 v)(fi 2 + i> 2 ) + sinh 2 /i sin^ v 



2 ,.12 



A; 2 (cosh ji cos v + 1) 



-—. — r- 7) n — : h F (tan<^, 

a(cosh /i + cos v) 2Md 2 sinh // sin v \ cosh + cos v 4 



(^sinh^sinh/sin^sin^)- 1 / 2 / dE x ^\^)$> x<j> {<\>") f 



dtj 
(310) 



H(s")=»" 



u(s")=u" 

P/i(s) / T>u(s) exp 



m(o)=m' 



1/(0)=!/' 



M 



(fi 2 + z> 2 ) + ^(cosh 7 /i - cos 2 v) 



+ad(cosh M -cos,)-— ^ g . nh ^ 
Unfortunately we cannot solve this path integral. 



h 2 I 'A 2 — \ + k\ cosh /x X 2 — j + k 2 cos z/ 



+ 



sin 2 v 



ds 
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We consider the potential (p = \J x 2 + y 2 



V 8 (x) 



a 



(312) 



This potential is the same potential as in 3.1.4 with an additional z-dependence via F(z). The 
corresponding path integral formulations and solutions in mutually orthogonal circular parabolic 
coordinates are therefore a straightforward combination of the results of 3.1.4 and Kp(z" , z';T). 



4. Summary. 

In the present work we have made the first step in the classification of all the presently possible 
path integral solutions of the Smorodinsky-Winternitz potentials in two and three dimensions. 
This approach as compared with the Schrodinger approach has advantages as well as disadvan- 
tages. The advantages are that the path integral gives a global view of the problem in question, 
whereas the Schrodinger equation only a local one. This property enables one to incorporate 
into the Feynman path integral many important features like topological effects, perturbation 
theory as well as the investigation of non-perturbative effects, a semi-classical expansion of the 
propagator, and many more. It is quite likely that the Feynman path integral is possibly the 
only consistent way to incorporate renormalization theory in a unified field theory. Therefore it 
is always important to investigate solvable quantum systems "from the point of view of fluctu- 
ating paths" . A further advantage is that the explicit computation of the propagator (the 
Feynman kernel), respectively the (energy-dependent) Green's function, gives simultaneously 
the spectral expansion into the wave-functions and the energy spectrum. 

A disadvantage is the fact that an explicit path integration can be done only in a limited 
number of coordinate systems, namely in polar, spherical, and some problems in parabolic 
coordinates, respectively coordinate systems which are related to them. Path integral evaluations 
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in, say, elliptic, paraboloidal and other parametric coordinate systems are not accessible. One 
may say that one of the very advantages of the Feynman path integral, i.e. its global approach 
to a quantum mechanical problem, leads in these cases to a disadvantage, because the actual 
path integration would require "addition theorems" in terms of special functions of parametric 
coordinate systems, i.e. their corresponding spherical functions of the Laplace operator. Such 
"addition theorems" are not to our disposal presently. However, it seems reasonable that a path 
integral identity in conical coordinates can be derived via the path integral on the sphere , 
where for the spherical functions a representation in terms of Jacobi functions is chosen. 

Explicitly solvable models always involve solutions in terms of hypergeometric and confluent 
hypergeometric functions. We have seen that in at least one coordinate system each of the dis- 
cussed Smorodinsky-Winternitz potentials can be solved in terms of these higher transcendental 
functions. In the parameter dependent coordinates the solution cannot be written in closed form. 
Therefore an investigation of the Feynman path integral in such coordinate systems is desirable 
and further studies along these lines will be subject to future work. This is not idle doing, then, 
for instance, the spheroidal coordinate system is only one into which the two-center problem 
separates. In future investigations we will try to discuss these more complicated problems. 

The present work does not claim a full analysis of the quantum mechanical Smorodinsky- 
Winternitz system. Such an analysis must include not only a comprehensive classification but 
also the determination of the interbasis coefficients of the wave-functions of the same energy. 
This requires the knowledge of the dynamical symmetry group which is responsible for the 
accidental degeneracy of the energy spectrum of a Smorodinsky-Winternitz potential. 

In due time we will continue our work along the lines described above, and we will also 
try to include an analysis in which quantum mechanical problems in spaces of (positive and 
negative) constant curvature are taken into account. This will include a path integral approach 
of parametric coordinate systems on the sphere and pseudo-sphere, the study of the Coulomb- 
and the Higgs oscillator problem in particular, and of Smorodinsky-Winternitz potentials in 
spaces of constant curvature in general. 
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Appendix. 

In this appendix we want to give a short discussion of the path integral problem of the anhar- 
monic sextic oscillator as encountered in parabolic coordinates for the Holt potential. Let us 
consider the one-dimensional path integral 



x(t")=x" 



K(x",x';T) = J Vx(t) exp J ^ 



x(t')=x' 



dt 



(313) 



Performing a combined space-time transformation with the new coordinate z = x 4 , i.e. x = 
F(z) = z 1 / 4 , and the new "pseudo-time" defined by Sj — Sj_i = 6j = 16x|x|_ 1 (tj — on the 
lattice, respectively s" = 16 jt x 6 (s)ds, we obtain the transformtion formulas (c.f. p8| , |40| , [59| , 
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K(x\x';T) = f°° ^ TE / n G(z",z';E) 



(314) 
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G(z",z';E) 



F'(z")F'(z') 



1/2 i 

h 



K(z",z';s")ds" , 



(315) 



and the kernel K(z", z'\ s") is given by [A 2 



(/3 2 + 15/4)16] 



z(s")- 



K(z",z',s") 



Vz(s) exp < — 



*(o)= 



M 



z 2 + 



2 1,2 



h z k 



+ 



E 



h 2 A 2 



32Mz I6z 3 / 2 2M 



ds 



(316) 



We see that this path integral cannot be evaluated for E ^ 0. This corresponds in the case of 
the Holt potential (and all related ones as well) that we can only "solve" it for the parabolic 
separation parameter equal to zero. 

If we now set E = the path integral K(s") is actually a Coulomb potential problem. For the 
Coulomb potential the Feynman kernel is not known in closed form, however the corresponding 
Green's function G^ C \E) is. Insertion in ( |315| ) gives the the zero-energy Green's function of 

m 



G(z",z';0) = F'(z")F'(z') 1/2 



z ,z; 



32 



(z'zT 



. 4 r(i + A + hk 2 /8Mio) 
Uj r(l + 2A) 

vff f MuJ 



M_ 



hk 2 /SMuj,\ 



(Muj 



(317) 



(318) 



Note that this formula is a generalization of the result of Steiner [119] where only one power 
potential term was taken into account. It is obvious that in a similar way an x 14 potential term 
can be considered instead of a x 6 , leading to a (zero-energy) Green's function of a radial harmonic 
oscillator. Furthermore, higher power potentials can be also treated by choosing another power 
dependence in the transformation x = F(z). 

Note that the zero-energy Green's function can also be evaluated in cases, where coordinate 
systems in flat space separate the Laplace equation (and not the Helmhotz, respectively the 
Schrodinger equation), i.e. coordinate systems which are R-separable, as for instance bispherical 
or toroidal coordinates. 
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